{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 79,
   "id": "377223cf",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "pd.options.display.float_format = '{:,.0f}'.format"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ef486e3a",
   "metadata": {},
   "source": [
    "## 버스승하차 데이터, 미세먼지 농도 데이터 불러오기"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "9fe0d1fa",
   "metadata": {},
   "outputs": [],
   "source": [
    "dust_df_GJ = pd.read_csv('dust_df_GJ.csv')\n",
    "bususer_GJ = pd.read_csv('bususer_GwangJin.csv')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 70,
   "id": "9aa77716",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>사용년월</th>\n",
       "      <th>표준버스정류장ID</th>\n",
       "      <th>버스정류장ARS번호</th>\n",
       "      <th>0</th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "      <th>5</th>\n",
       "      <th>6</th>\n",
       "      <th>...</th>\n",
       "      <th>17</th>\n",
       "      <th>18</th>\n",
       "      <th>19</th>\n",
       "      <th>20</th>\n",
       "      <th>21</th>\n",
       "      <th>22</th>\n",
       "      <th>23</th>\n",
       "      <th>행정동</th>\n",
       "      <th>역명</th>\n",
       "      <th>시간합</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2914</th>\n",
       "      <td>11</td>\n",
       "      <td>104900098</td>\n",
       "      <td>5710</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>48</td>\n",
       "      <td>...</td>\n",
       "      <td>161</td>\n",
       "      <td>177</td>\n",
       "      <td>116</td>\n",
       "      <td>95</td>\n",
       "      <td>76</td>\n",
       "      <td>40</td>\n",
       "      <td>12</td>\n",
       "      <td>구의1동</td>\n",
       "      <td>미니스톱구의중앙점앞</td>\n",
       "      <td>3149</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>423</th>\n",
       "      <td>4</td>\n",
       "      <td>104000037</td>\n",
       "      <td>5130</td>\n",
       "      <td>326</td>\n",
       "      <td>5</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>37</td>\n",
       "      <td>195</td>\n",
       "      <td>...</td>\n",
       "      <td>617</td>\n",
       "      <td>704</td>\n",
       "      <td>666</td>\n",
       "      <td>480</td>\n",
       "      <td>487</td>\n",
       "      <td>466</td>\n",
       "      <td>339</td>\n",
       "      <td>중곡3동</td>\n",
       "      <td>면곡시장앞</td>\n",
       "      <td>8244</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>397</th>\n",
       "      <td>2</td>\n",
       "      <td>104000035</td>\n",
       "      <td>5128</td>\n",
       "      <td>2</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>337</td>\n",
       "      <td>592</td>\n",
       "      <td>907</td>\n",
       "      <td>...</td>\n",
       "      <td>1030</td>\n",
       "      <td>1204</td>\n",
       "      <td>877</td>\n",
       "      <td>669</td>\n",
       "      <td>731</td>\n",
       "      <td>518</td>\n",
       "      <td>142</td>\n",
       "      <td>중곡3동</td>\n",
       "      <td>중곡3동주민센터앞</td>\n",
       "      <td>18573</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>937</th>\n",
       "      <td>2</td>\n",
       "      <td>104000095</td>\n",
       "      <td>5188</td>\n",
       "      <td>416</td>\n",
       "      <td>86</td>\n",
       "      <td>13</td>\n",
       "      <td>0</td>\n",
       "      <td>296</td>\n",
       "      <td>705</td>\n",
       "      <td>1743</td>\n",
       "      <td>...</td>\n",
       "      <td>3306</td>\n",
       "      <td>3619</td>\n",
       "      <td>2596</td>\n",
       "      <td>1881</td>\n",
       "      <td>2047</td>\n",
       "      <td>1593</td>\n",
       "      <td>899</td>\n",
       "      <td>자양1동</td>\n",
       "      <td>자양동잠실대교북단</td>\n",
       "      <td>50839</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1648</th>\n",
       "      <td>5</td>\n",
       "      <td>104000167</td>\n",
       "      <td>5260</td>\n",
       "      <td>123</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>370</td>\n",
       "      <td>909</td>\n",
       "      <td>1372</td>\n",
       "      <td>...</td>\n",
       "      <td>1856</td>\n",
       "      <td>1858</td>\n",
       "      <td>1312</td>\n",
       "      <td>1017</td>\n",
       "      <td>867</td>\n",
       "      <td>673</td>\n",
       "      <td>302</td>\n",
       "      <td>중곡2동</td>\n",
       "      <td>중곡동입구</td>\n",
       "      <td>32811</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5 rows × 30 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      사용년월  표준버스정류장ID  버스정류장ARS번호    0   1   2  3    4    5     6  ...    17  \\\n",
       "2914    11  104900098        5710    0   0   0  0    0    0    48  ...   161   \n",
       "423      4  104000037        5130  326   5   0  0    0   37   195  ...   617   \n",
       "397      2  104000035        5128    2   0   0  0  337  592   907  ...  1030   \n",
       "937      2  104000095        5188  416  86  13  0  296  705  1743  ...  3306   \n",
       "1648     5  104000167        5260  123   0   0  0  370  909  1372  ...  1856   \n",
       "\n",
       "        18    19    20    21    22   23   행정동          역명    시간합  \n",
       "2914   177   116    95    76    40   12  구의1동  미니스톱구의중앙점앞   3149  \n",
       "423    704   666   480   487   466  339  중곡3동       면곡시장앞   8244  \n",
       "397   1204   877   669   731   518  142  중곡3동   중곡3동주민센터앞  18573  \n",
       "937   3619  2596  1881  2047  1593  899  자양1동   자양동잠실대교북단  50839  \n",
       "1648  1858  1312  1017   867   673  302  중곡2동       중곡동입구  32811  \n",
       "\n",
       "[5 rows x 30 columns]"
      ]
     },
     "execution_count": 70,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bususer_GJ.sample(5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "28aeac4b",
   "metadata": {},
   "outputs": [],
   "source": [
    "dust_df_GJ = dust_df_GJ.iloc[:,3:]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "30e96183",
   "metadata": {},
   "outputs": [],
   "source": [
    "bususer_GJ = bususer_GJ.iloc[:,1:]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 72,
   "id": "20f066c2",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>초미세먼지 보정(㎍/㎥)</th>\n",
       "      <th>미세먼지 보정(㎍/㎥)</th>\n",
       "      <th>date</th>\n",
       "      <th>month</th>\n",
       "      <th>day</th>\n",
       "      <th>hour</th>\n",
       "      <th>행정동</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>31537</th>\n",
       "      <td>28.00</td>\n",
       "      <td>40.00</td>\n",
       "      <td>2022-04-17 14:00:00</td>\n",
       "      <td>4</td>\n",
       "      <td>17</td>\n",
       "      <td>14</td>\n",
       "      <td>Guui3-dong</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74228</th>\n",
       "      <td>4.00</td>\n",
       "      <td>8.00</td>\n",
       "      <td>2022-09-18 23:00:00</td>\n",
       "      <td>9</td>\n",
       "      <td>18</td>\n",
       "      <td>23</td>\n",
       "      <td>Gwangjang-dong</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>100142</th>\n",
       "      <td>27.00</td>\n",
       "      <td>29.00</td>\n",
       "      <td>2022-12-25 13:00:00</td>\n",
       "      <td>12</td>\n",
       "      <td>25</td>\n",
       "      <td>13</td>\n",
       "      <td>Guui3-dong</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>58383</th>\n",
       "      <td>9.00</td>\n",
       "      <td>7.00</td>\n",
       "      <td>2022-07-27 11:00:00</td>\n",
       "      <td>7</td>\n",
       "      <td>27</td>\n",
       "      <td>11</td>\n",
       "      <td>Neung-dong</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>90315</th>\n",
       "      <td>15.00</td>\n",
       "      <td>27.00</td>\n",
       "      <td>2022-11-07 13:00:00</td>\n",
       "      <td>11</td>\n",
       "      <td>7</td>\n",
       "      <td>13</td>\n",
       "      <td>Gwangjang-dong</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        초미세먼지 보정(㎍/㎥)  미세먼지 보정(㎍/㎥)                 date  month  day  hour  \\\n",
       "31537           28.00         40.00  2022-04-17 14:00:00      4   17    14   \n",
       "74228            4.00          8.00  2022-09-18 23:00:00      9   18    23   \n",
       "100142          27.00         29.00  2022-12-25 13:00:00     12   25    13   \n",
       "58383            9.00          7.00  2022-07-27 11:00:00      7   27    11   \n",
       "90315           15.00         27.00  2022-11-07 13:00:00     11    7    13   \n",
       "\n",
       "                   행정동  \n",
       "31537       Guui3-dong  \n",
       "74228   Gwangjang-dong  \n",
       "100142      Guui3-dong  \n",
       "58383       Neung-dong  \n",
       "90315   Gwangjang-dong  "
      ]
     },
     "execution_count": 72,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dust_df_GJ.iloc[:,3:].sample(5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "id": "b6da0671",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "19.669715979254136"
      ]
     },
     "execution_count": 36,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dust_df_GJ['미세먼지 보정(㎍/㎥)'].mean()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "379914f4",
   "metadata": {},
   "outputs": [],
   "source": [
    "#사용년월 형식 변경\n",
    "bususer_GJ['사용년월'] = bususer_GJ['사용년월']%202200"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "1555bfbc",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>사용년월</th>\n",
       "      <th>표준버스정류장ID</th>\n",
       "      <th>버스정류장ARS번호</th>\n",
       "      <th>0</th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "      <th>5</th>\n",
       "      <th>6</th>\n",
       "      <th>...</th>\n",
       "      <th>16</th>\n",
       "      <th>17</th>\n",
       "      <th>18</th>\n",
       "      <th>19</th>\n",
       "      <th>20</th>\n",
       "      <th>21</th>\n",
       "      <th>22</th>\n",
       "      <th>23</th>\n",
       "      <th>행정동</th>\n",
       "      <th>역명</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>101000252</td>\n",
       "      <td>5118</td>\n",
       "      <td>140</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>418</td>\n",
       "      <td>966</td>\n",
       "      <td>1666</td>\n",
       "      <td>...</td>\n",
       "      <td>2493</td>\n",
       "      <td>2920</td>\n",
       "      <td>2617</td>\n",
       "      <td>2103</td>\n",
       "      <td>1503</td>\n",
       "      <td>1514</td>\n",
       "      <td>713</td>\n",
       "      <td>363</td>\n",
       "      <td>중곡2동</td>\n",
       "      <td>신성시장앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>101000252</td>\n",
       "      <td>5118</td>\n",
       "      <td>81</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>380</td>\n",
       "      <td>929</td>\n",
       "      <td>1520</td>\n",
       "      <td>...</td>\n",
       "      <td>2302</td>\n",
       "      <td>2415</td>\n",
       "      <td>2501</td>\n",
       "      <td>1937</td>\n",
       "      <td>1307</td>\n",
       "      <td>1197</td>\n",
       "      <td>729</td>\n",
       "      <td>271</td>\n",
       "      <td>중곡2동</td>\n",
       "      <td>신성시장앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>101000252</td>\n",
       "      <td>5118</td>\n",
       "      <td>163</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>481</td>\n",
       "      <td>1017</td>\n",
       "      <td>1944</td>\n",
       "      <td>...</td>\n",
       "      <td>2594</td>\n",
       "      <td>2782</td>\n",
       "      <td>2823</td>\n",
       "      <td>2232</td>\n",
       "      <td>1504</td>\n",
       "      <td>1164</td>\n",
       "      <td>863</td>\n",
       "      <td>406</td>\n",
       "      <td>중곡2동</td>\n",
       "      <td>신성시장앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>4</td>\n",
       "      <td>101000252</td>\n",
       "      <td>5118</td>\n",
       "      <td>220</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>459</td>\n",
       "      <td>989</td>\n",
       "      <td>2099</td>\n",
       "      <td>...</td>\n",
       "      <td>2709</td>\n",
       "      <td>2942</td>\n",
       "      <td>2868</td>\n",
       "      <td>2358</td>\n",
       "      <td>1613</td>\n",
       "      <td>1259</td>\n",
       "      <td>891</td>\n",
       "      <td>408</td>\n",
       "      <td>중곡2동</td>\n",
       "      <td>신성시장앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>5</td>\n",
       "      <td>101000252</td>\n",
       "      <td>5118</td>\n",
       "      <td>221</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>477</td>\n",
       "      <td>1051</td>\n",
       "      <td>2200</td>\n",
       "      <td>...</td>\n",
       "      <td>2877</td>\n",
       "      <td>3103</td>\n",
       "      <td>3261</td>\n",
       "      <td>2327</td>\n",
       "      <td>1859</td>\n",
       "      <td>1351</td>\n",
       "      <td>930</td>\n",
       "      <td>413</td>\n",
       "      <td>중곡2동</td>\n",
       "      <td>신성시장앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3223</th>\n",
       "      <td>8</td>\n",
       "      <td>104900129</td>\n",
       "      <td>5999</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1710</td>\n",
       "      <td>...</td>\n",
       "      <td>3057</td>\n",
       "      <td>4067</td>\n",
       "      <td>5479</td>\n",
       "      <td>3945</td>\n",
       "      <td>2880</td>\n",
       "      <td>2198</td>\n",
       "      <td>2003</td>\n",
       "      <td>787</td>\n",
       "      <td>구의3동</td>\n",
       "      <td>강변역.테크노마트앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3224</th>\n",
       "      <td>9</td>\n",
       "      <td>104900129</td>\n",
       "      <td>5999</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1768</td>\n",
       "      <td>...</td>\n",
       "      <td>3067</td>\n",
       "      <td>4046</td>\n",
       "      <td>5505</td>\n",
       "      <td>3756</td>\n",
       "      <td>2931</td>\n",
       "      <td>2314</td>\n",
       "      <td>2121</td>\n",
       "      <td>721</td>\n",
       "      <td>구의3동</td>\n",
       "      <td>강변역.테크노마트앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3225</th>\n",
       "      <td>10</td>\n",
       "      <td>104900129</td>\n",
       "      <td>5999</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1673</td>\n",
       "      <td>...</td>\n",
       "      <td>3191</td>\n",
       "      <td>4049</td>\n",
       "      <td>5726</td>\n",
       "      <td>4008</td>\n",
       "      <td>3136</td>\n",
       "      <td>2346</td>\n",
       "      <td>2175</td>\n",
       "      <td>782</td>\n",
       "      <td>구의3동</td>\n",
       "      <td>강변역.테크노마트앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3226</th>\n",
       "      <td>11</td>\n",
       "      <td>104900129</td>\n",
       "      <td>5999</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1634</td>\n",
       "      <td>...</td>\n",
       "      <td>3112</td>\n",
       "      <td>4182</td>\n",
       "      <td>5907</td>\n",
       "      <td>4092</td>\n",
       "      <td>2875</td>\n",
       "      <td>2363</td>\n",
       "      <td>1981</td>\n",
       "      <td>756</td>\n",
       "      <td>구의3동</td>\n",
       "      <td>강변역.테크노마트앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3227</th>\n",
       "      <td>12</td>\n",
       "      <td>104900129</td>\n",
       "      <td>5999</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1619</td>\n",
       "      <td>...</td>\n",
       "      <td>3410</td>\n",
       "      <td>4264</td>\n",
       "      <td>5620</td>\n",
       "      <td>3747</td>\n",
       "      <td>2827</td>\n",
       "      <td>2229</td>\n",
       "      <td>2065</td>\n",
       "      <td>777</td>\n",
       "      <td>구의3동</td>\n",
       "      <td>강변역.테크노마트앞</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3228 rows × 29 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      사용년월  표준버스정류장ID  버스정류장ARS번호    0  1  2  3    4     5     6  ...    16  \\\n",
       "0        1  101000252        5118  140  0  0  0  418   966  1666  ...  2493   \n",
       "1        2  101000252        5118   81  0  0  0  380   929  1520  ...  2302   \n",
       "2        3  101000252        5118  163  0  0  0  481  1017  1944  ...  2594   \n",
       "3        4  101000252        5118  220  0  0  0  459   989  2099  ...  2709   \n",
       "4        5  101000252        5118  221  0  0  0  477  1051  2200  ...  2877   \n",
       "...    ...        ...         ...  ... .. .. ..  ...   ...   ...  ...   ...   \n",
       "3223     8  104900129        5999    0  0  0  0    0     0  1710  ...  3057   \n",
       "3224     9  104900129        5999    0  0  0  0    0     0  1768  ...  3067   \n",
       "3225    10  104900129        5999    0  0  0  0    0     0  1673  ...  3191   \n",
       "3226    11  104900129        5999    0  0  0  0    0     0  1634  ...  3112   \n",
       "3227    12  104900129        5999    0  0  0  0    0     0  1619  ...  3410   \n",
       "\n",
       "        17    18    19    20    21    22   23   행정동          역명  \n",
       "0     2920  2617  2103  1503  1514   713  363  중곡2동       신성시장앞  \n",
       "1     2415  2501  1937  1307  1197   729  271  중곡2동       신성시장앞  \n",
       "2     2782  2823  2232  1504  1164   863  406  중곡2동       신성시장앞  \n",
       "3     2942  2868  2358  1613  1259   891  408  중곡2동       신성시장앞  \n",
       "4     3103  3261  2327  1859  1351   930  413  중곡2동       신성시장앞  \n",
       "...    ...   ...   ...   ...   ...   ...  ...   ...         ...  \n",
       "3223  4067  5479  3945  2880  2198  2003  787  구의3동  강변역.테크노마트앞  \n",
       "3224  4046  5505  3756  2931  2314  2121  721  구의3동  강변역.테크노마트앞  \n",
       "3225  4049  5726  4008  3136  2346  2175  782  구의3동  강변역.테크노마트앞  \n",
       "3226  4182  5907  4092  2875  2363  1981  756  구의3동  강변역.테크노마트앞  \n",
       "3227  4264  5620  3747  2827  2229  2065  777  구의3동  강변역.테크노마트앞  \n",
       "\n",
       "[3228 rows x 29 columns]"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bususer_GJ"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "6da696e5",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "114d0eee",
   "metadata": {},
   "source": [
    "## groupby 데이터프레임 만들기"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "0101d1f1",
   "metadata": {},
   "outputs": [],
   "source": [
    "# 자양2동은 오류데이터가 많으므로 지리적으로 가까운 자양3동의 데이터로 대체\n",
    "guui1 = dust_df_GJ[dust_df_GJ['행정동']=='Guui1-dong'][['초미세먼지 보정(㎍/㎥)', '미세먼지 보정(㎍/㎥)', 'month', 'hour']].groupby(['month', 'hour']).mean().reset_index()\n",
    "guui2 = dust_df_GJ[dust_df_GJ['행정동']=='Guui2-dong'][['초미세먼지 보정(㎍/㎥)', '미세먼지 보정(㎍/㎥)', 'month', 'hour']].groupby(['month', 'hour']).mean().reset_index()\n",
    "guui3 = dust_df_GJ[dust_df_GJ['행정동']=='Guui3-dong'][['초미세먼지 보정(㎍/㎥)', '미세먼지 보정(㎍/㎥)', 'month', 'hour']].groupby(['month', 'hour']).mean().reset_index()\n",
    "jayang1 = dust_df_GJ[dust_df_GJ['행정동']=='Jayang1-dong'][['초미세먼지 보정(㎍/㎥)', '미세먼지 보정(㎍/㎥)', 'month', 'hour']].groupby(['month', 'hour']).mean().reset_index()\n",
    "jayang2 = dust_df_GJ[dust_df_GJ['행정동']=='Jayang3-dong'][['초미세먼지 보정(㎍/㎥)', '미세먼지 보정(㎍/㎥)', 'month', 'hour']].groupby(['month', 'hour']).mean().reset_index()\n",
    "jayang3 = dust_df_GJ[dust_df_GJ['행정동']=='Jayang3-dong'][['초미세먼지 보정(㎍/㎥)', '미세먼지 보정(㎍/㎥)', 'month', 'hour']].groupby(['month', 'hour']).mean().reset_index()\n",
    "jayang4 = dust_df_GJ[dust_df_GJ['행정동']=='Jayang4-dong'][['초미세먼지 보정(㎍/㎥)', '미세먼지 보정(㎍/㎥)', 'month', 'hour']].groupby(['month', 'hour']).mean().reset_index()\n",
    "junggok1 = dust_df_GJ[dust_df_GJ['행정동']=='Junggok1-dong'][['초미세먼지 보정(㎍/㎥)', '미세먼지 보정(㎍/㎥)', 'month', 'hour']].groupby(['month', 'hour']).mean().reset_index()\n",
    "junggok2 = dust_df_GJ[dust_df_GJ['행정동']=='Junggok2-dong'][['초미세먼지 보정(㎍/㎥)', '미세먼지 보정(㎍/㎥)', 'month', 'hour']].groupby(['month', 'hour']).mean().reset_index()\n",
    "junggok3 = dust_df_GJ[dust_df_GJ['행정동']=='Junggok3-dong'][['초미세먼지 보정(㎍/㎥)', '미세먼지 보정(㎍/㎥)', 'month', 'hour']].groupby(['month', 'hour']).mean().reset_index()\n",
    "junggok4 = dust_df_GJ[dust_df_GJ['행정동']=='Junggok4-dong'][['초미세먼지 보정(㎍/㎥)', '미세먼지 보정(㎍/㎥)', 'month', 'hour']].groupby(['month', 'hour']).mean().reset_index()\n",
    "neung = dust_df_GJ[dust_df_GJ['행정동']=='Neung-dong'][['초미세먼지 보정(㎍/㎥)', '미세먼지 보정(㎍/㎥)', 'month', 'hour']].groupby(['month', 'hour']).mean().reset_index()\n",
    "gunja = dust_df_GJ[dust_df_GJ['행정동']=='Gunja-dong'][['초미세먼지 보정(㎍/㎥)', '미세먼지 보정(㎍/㎥)', 'month', 'hour']].groupby(['month', 'hour']).mean().reset_index()\n",
    "hwayang = dust_df_GJ[dust_df_GJ['행정동']=='Hwayang-dong'][['초미세먼지 보정(㎍/㎥)', '미세먼지 보정(㎍/㎥)', 'month', 'hour']].groupby(['month', 'hour']).mean().reset_index()\n",
    "gwangjang = dust_df_GJ[dust_df_GJ['행정동']=='Gwangjang-dong'][['초미세먼지 보정(㎍/㎥)', '미세먼지 보정(㎍/㎥)', 'month', 'hour']].groupby(['month', 'hour']).mean().reset_index()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "3d55ffc5",
   "metadata": {},
   "outputs": [],
   "source": [
    "dong_list = []\n",
    "dong_list.extend([guui1, guui2, guui3, jayang1, jayang2, jayang3, jayang4, junggok1, junggok2, junggok3, junggok4, neung, gunja, hwayang, gwangjang])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "f15f7588",
   "metadata": {},
   "outputs": [],
   "source": [
    "dong_dict = {'구의1동':0, '구의2동':1, '구의3동':2, '자양1동':3, '자양2동':4, '자양3동':5, \n",
    "             '자양4동':6, '중곡1동':7, '중곡2동':8, '중곡3동':9, '중곡4동':10, '능동':11, \n",
    "             '군자동':12, '화양동':13, '광장동':14}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "f5fb6d0e",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>month</th>\n",
       "      <th>hour</th>\n",
       "      <th>초미세먼지 보정(㎍/㎥)</th>\n",
       "      <th>미세먼지 보정(㎍/㎥)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>264</th>\n",
       "      <td>12</td>\n",
       "      <td>0</td>\n",
       "      <td>15.400000</td>\n",
       "      <td>17.900000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>265</th>\n",
       "      <td>12</td>\n",
       "      <td>1</td>\n",
       "      <td>7.000000</td>\n",
       "      <td>9.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>266</th>\n",
       "      <td>12</td>\n",
       "      <td>2</td>\n",
       "      <td>6.000000</td>\n",
       "      <td>6.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>267</th>\n",
       "      <td>12</td>\n",
       "      <td>4</td>\n",
       "      <td>14.888889</td>\n",
       "      <td>18.222222</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>268</th>\n",
       "      <td>12</td>\n",
       "      <td>5</td>\n",
       "      <td>16.900000</td>\n",
       "      <td>20.200000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>269</th>\n",
       "      <td>12</td>\n",
       "      <td>6</td>\n",
       "      <td>18.800000</td>\n",
       "      <td>23.800000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>270</th>\n",
       "      <td>12</td>\n",
       "      <td>7</td>\n",
       "      <td>15.500000</td>\n",
       "      <td>17.500000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>271</th>\n",
       "      <td>12</td>\n",
       "      <td>8</td>\n",
       "      <td>17.400000</td>\n",
       "      <td>20.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>272</th>\n",
       "      <td>12</td>\n",
       "      <td>9</td>\n",
       "      <td>16.500000</td>\n",
       "      <td>20.100000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>273</th>\n",
       "      <td>12</td>\n",
       "      <td>10</td>\n",
       "      <td>18.100000</td>\n",
       "      <td>21.700000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>274</th>\n",
       "      <td>12</td>\n",
       "      <td>11</td>\n",
       "      <td>15.800000</td>\n",
       "      <td>18.700000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>275</th>\n",
       "      <td>12</td>\n",
       "      <td>12</td>\n",
       "      <td>15.600000</td>\n",
       "      <td>18.500000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>276</th>\n",
       "      <td>12</td>\n",
       "      <td>13</td>\n",
       "      <td>16.300000</td>\n",
       "      <td>18.600000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>277</th>\n",
       "      <td>12</td>\n",
       "      <td>14</td>\n",
       "      <td>16.400000</td>\n",
       "      <td>18.500000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>278</th>\n",
       "      <td>12</td>\n",
       "      <td>15</td>\n",
       "      <td>15.700000</td>\n",
       "      <td>17.500000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>279</th>\n",
       "      <td>12</td>\n",
       "      <td>16</td>\n",
       "      <td>13.666667</td>\n",
       "      <td>15.666667</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>280</th>\n",
       "      <td>12</td>\n",
       "      <td>17</td>\n",
       "      <td>17.600000</td>\n",
       "      <td>20.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>281</th>\n",
       "      <td>12</td>\n",
       "      <td>18</td>\n",
       "      <td>17.800000</td>\n",
       "      <td>19.400000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>282</th>\n",
       "      <td>12</td>\n",
       "      <td>19</td>\n",
       "      <td>17.300000</td>\n",
       "      <td>20.500000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>283</th>\n",
       "      <td>12</td>\n",
       "      <td>20</td>\n",
       "      <td>15.300000</td>\n",
       "      <td>17.200000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>284</th>\n",
       "      <td>12</td>\n",
       "      <td>21</td>\n",
       "      <td>16.700000</td>\n",
       "      <td>19.200000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>285</th>\n",
       "      <td>12</td>\n",
       "      <td>22</td>\n",
       "      <td>18.400000</td>\n",
       "      <td>22.200000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>286</th>\n",
       "      <td>12</td>\n",
       "      <td>23</td>\n",
       "      <td>14.900000</td>\n",
       "      <td>17.200000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     month  hour  초미세먼지 보정(㎍/㎥)  미세먼지 보정(㎍/㎥)\n",
       "264     12     0      15.400000     17.900000\n",
       "265     12     1       7.000000      9.000000\n",
       "266     12     2       6.000000      6.000000\n",
       "267     12     4      14.888889     18.222222\n",
       "268     12     5      16.900000     20.200000\n",
       "269     12     6      18.800000     23.800000\n",
       "270     12     7      15.500000     17.500000\n",
       "271     12     8      17.400000     20.000000\n",
       "272     12     9      16.500000     20.100000\n",
       "273     12    10      18.100000     21.700000\n",
       "274     12    11      15.800000     18.700000\n",
       "275     12    12      15.600000     18.500000\n",
       "276     12    13      16.300000     18.600000\n",
       "277     12    14      16.400000     18.500000\n",
       "278     12    15      15.700000     17.500000\n",
       "279     12    16      13.666667     15.666667\n",
       "280     12    17      17.600000     20.000000\n",
       "281     12    18      17.800000     19.400000\n",
       "282     12    19      17.300000     20.500000\n",
       "283     12    20      15.300000     17.200000\n",
       "284     12    21      16.700000     19.200000\n",
       "285     12    22      18.400000     22.200000\n",
       "286     12    23      14.900000     17.200000"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "junggok1[junggok1['month'] == 12]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0d67327f",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0e5f35dc",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "3ef260f8",
   "metadata": {},
   "source": [
    "## 버스승하차 수와 미세먼지 농도 곱하기"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "7d759800",
   "metadata": {},
   "outputs": [],
   "source": [
    "bus_dust10 = bususer_GJ.copy()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "fa9acd6c",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0\n",
      "1\n",
      "2\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "중곡1동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "3\n",
      "4\n",
      "5\n",
      "6\n",
      "7\n",
      "8\n",
      "9\n",
      "10\n",
      "11\n",
      "12\n",
      "13\n",
      "14\n",
      "15\n",
      "16\n",
      "17\n",
      "18\n",
      "19\n",
      "20\n",
      "21\n",
      "22\n",
      "23\n"
     ]
    }
   ],
   "source": [
    "for j in range(24):\n",
    "    hour = []\n",
    "    for i, dong in enumerate(bus_dust10['행정동']):\n",
    "        temp_bus = bususer_GJ.iloc[:,[0,3+j]]\n",
    "        temp_dust = dong_list[dong_dict[dong]]\n",
    "        bus = temp_bus.iloc[i,1]\n",
    "        temp1 = temp_dust[temp_dust['month'] == temp_bus.iloc[i,0]]\n",
    "        temp2 = temp1[temp1['hour'] == j]\n",
    "        \n",
    "        if temp2.shape[0] == 0:\n",
    "            dust=0\n",
    "            print(dong)\n",
    "            print(temp_bus.iloc[i,0])\n",
    "        else:\n",
    "            dust=temp2.iloc[0,3]\n",
    "        hour.append(bus*dust)\n",
    "    bus_dust10[str(j)] = hour\n",
    "    print(j)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 80,
   "id": "7b15ff8a",
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>사용년월</th>\n",
       "      <th>표준버스정류장ID</th>\n",
       "      <th>버스정류장ARS번호</th>\n",
       "      <th>0</th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "      <th>5</th>\n",
       "      <th>6</th>\n",
       "      <th>...</th>\n",
       "      <th>17</th>\n",
       "      <th>18</th>\n",
       "      <th>19</th>\n",
       "      <th>20</th>\n",
       "      <th>21</th>\n",
       "      <th>22</th>\n",
       "      <th>23</th>\n",
       "      <th>행정동</th>\n",
       "      <th>역명</th>\n",
       "      <th>시간합</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1886</th>\n",
       "      <td>3</td>\n",
       "      <td>104000299</td>\n",
       "      <td>5008</td>\n",
       "      <td>326</td>\n",
       "      <td>44</td>\n",
       "      <td>20</td>\n",
       "      <td>46</td>\n",
       "      <td>206</td>\n",
       "      <td>900</td>\n",
       "      <td>886</td>\n",
       "      <td>...</td>\n",
       "      <td>3,001</td>\n",
       "      <td>3,528</td>\n",
       "      <td>2,278</td>\n",
       "      <td>1,509</td>\n",
       "      <td>1,170</td>\n",
       "      <td>1,039</td>\n",
       "      <td>721</td>\n",
       "      <td>능동</td>\n",
       "      <td>중곡동입구삼거리</td>\n",
       "      <td>704,098</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2150</th>\n",
       "      <td>3</td>\n",
       "      <td>104900020</td>\n",
       "      <td>5578</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>95</td>\n",
       "      <td>...</td>\n",
       "      <td>593</td>\n",
       "      <td>496</td>\n",
       "      <td>638</td>\n",
       "      <td>299</td>\n",
       "      <td>338</td>\n",
       "      <td>333</td>\n",
       "      <td>108</td>\n",
       "      <td>광장동</td>\n",
       "      <td>광장중학교앞</td>\n",
       "      <td>138,528</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2267</th>\n",
       "      <td>12</td>\n",
       "      <td>104900029</td>\n",
       "      <td>5559</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>20</td>\n",
       "      <td>...</td>\n",
       "      <td>219</td>\n",
       "      <td>254</td>\n",
       "      <td>271</td>\n",
       "      <td>130</td>\n",
       "      <td>123</td>\n",
       "      <td>82</td>\n",
       "      <td>0</td>\n",
       "      <td>중곡4동</td>\n",
       "      <td>긴고랑입구</td>\n",
       "      <td>43,106</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>724</th>\n",
       "      <td>5</td>\n",
       "      <td>104000062</td>\n",
       "      <td>5155</td>\n",
       "      <td>25</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>36</td>\n",
       "      <td>152</td>\n",
       "      <td>228</td>\n",
       "      <td>...</td>\n",
       "      <td>351</td>\n",
       "      <td>420</td>\n",
       "      <td>297</td>\n",
       "      <td>249</td>\n",
       "      <td>241</td>\n",
       "      <td>267</td>\n",
       "      <td>167</td>\n",
       "      <td>구의2동</td>\n",
       "      <td>래미안구의파크스위트</td>\n",
       "      <td>223,185</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1178</th>\n",
       "      <td>3</td>\n",
       "      <td>104000123</td>\n",
       "      <td>5216</td>\n",
       "      <td>412</td>\n",
       "      <td>94</td>\n",
       "      <td>43</td>\n",
       "      <td>43</td>\n",
       "      <td>95</td>\n",
       "      <td>425</td>\n",
       "      <td>937</td>\n",
       "      <td>...</td>\n",
       "      <td>3,472</td>\n",
       "      <td>4,090</td>\n",
       "      <td>2,660</td>\n",
       "      <td>2,056</td>\n",
       "      <td>2,302</td>\n",
       "      <td>2,099</td>\n",
       "      <td>1,021</td>\n",
       "      <td>자양4동</td>\n",
       "      <td>노룬산시장</td>\n",
       "      <td>715,153</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5 rows × 21 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      사용년월  표준버스정류장ID  버스정류장ARS번호   0  1  2  3   4   5   6  ...    17    18  \\\n",
       "1886     3  104000299        5008 326 44 20 46 206 900 886  ... 3,001 3,528   \n",
       "2150     3  104900020        5578   0  0  0  0   0   0  95  ...   593   496   \n",
       "2267    12  104900029        5559   0  0  0  0   0   0  20  ...   219   254   \n",
       "724      5  104000062        5155  25  0  0  0  36 152 228  ...   351   420   \n",
       "1178     3  104000123        5216 412 94 43 43  95 425 937  ... 3,472 4,090   \n",
       "\n",
       "        19    20    21    22    23   행정동          역명     시간합  \n",
       "1886 2,278 1,509 1,170 1,039   721    능동    중곡동입구삼거리 704,098  \n",
       "2150   638   299   338   333   108   광장동      광장중학교앞 138,528  \n",
       "2267   271   130   123    82     0  중곡4동       긴고랑입구  43,106  \n",
       "724    297   249   241   267   167  구의2동  래미안구의파크스위트 223,185  \n",
       "1178 2,660 2,056 2,302 2,099 1,021  자양4동       노룬산시장 715,153  \n",
       "\n",
       "[5 rows x 21 columns]"
      ]
     },
     "execution_count": 80,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bus_dust10_1 = bus_dust10.copy()\n",
    "bus_dust10_1.iloc[:,3:27]  = bus_dust10.iloc[:,3:27]/dust_df_GJ['미세먼지 보정(㎍/㎥)'].mean()\n",
    "bus_dust10_1.drop(['7','8','9','10','11', '12', '13', '14', '15'], axis=1).sample(5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "0527a069",
   "metadata": {},
   "outputs": [],
   "source": [
    "bus_dust10.to_csv('bus_dust10.csv')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "467348b1",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "2383fc14",
   "metadata": {},
   "source": [
    "## 버스승하차 수와 초미세먼지 농도 곱하기"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "0190d527",
   "metadata": {},
   "outputs": [],
   "source": [
    "bus_dust2_5 = bususer_GJ.copy()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "38e55fae",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0\n",
      "1\n",
      "2\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "중곡1동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡1동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "자양1동\n",
      "6\n",
      "자양1동\n",
      "6\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡4동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "중곡3동\n",
      "12\n",
      "3\n",
      "4\n",
      "5\n",
      "6\n",
      "7\n",
      "8\n",
      "9\n",
      "10\n",
      "11\n",
      "12\n",
      "13\n",
      "14\n",
      "15\n",
      "16\n",
      "17\n",
      "18\n",
      "19\n",
      "20\n",
      "21\n",
      "22\n",
      "23\n"
     ]
    }
   ],
   "source": [
    "#자양1동 -> 6월의 3시데이터 누락, 중곡1~4동 12월의 3시데이터 누락\n",
    "for j in range(24):\n",
    "    hour = []\n",
    "    for i, dong in enumerate(bus_dust2_5['행정동']):\n",
    "        temp_bus = bususer_GJ.iloc[:,[0,3+j]]\n",
    "        temp_dust = dong_list[dong_dict[dong]]\n",
    "        bus = temp_bus.iloc[i,1]\n",
    "        temp1 = temp_dust[temp_dust['month'] == temp_bus.iloc[i,0]]\n",
    "        temp2 = temp1[temp1['hour'] == j]\n",
    "        \n",
    "        if temp2.shape[0] == 0:\n",
    "            dust=0\n",
    "            print(dong)\n",
    "            print(temp_bus.iloc[i,0])\n",
    "        else:\n",
    "            dust=temp2.iloc[0,2]\n",
    "        hour.append(bus*dust)\n",
    "    bus_dust2_5[str(j)] = hour\n",
    "    print(j)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "2261578e",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>사용년월</th>\n",
       "      <th>표준버스정류장ID</th>\n",
       "      <th>버스정류장ARS번호</th>\n",
       "      <th>0</th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "      <th>5</th>\n",
       "      <th>6</th>\n",
       "      <th>...</th>\n",
       "      <th>16</th>\n",
       "      <th>17</th>\n",
       "      <th>18</th>\n",
       "      <th>19</th>\n",
       "      <th>20</th>\n",
       "      <th>21</th>\n",
       "      <th>22</th>\n",
       "      <th>23</th>\n",
       "      <th>행정동</th>\n",
       "      <th>역명</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>101000252</td>\n",
       "      <td>5118</td>\n",
       "      <td>3944.137931</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>11022.814815</td>\n",
       "      <td>28014.00000</td>\n",
       "      <td>49520.413793</td>\n",
       "      <td>...</td>\n",
       "      <td>73242.620690</td>\n",
       "      <td>81659.310345</td>\n",
       "      <td>76253.965517</td>\n",
       "      <td>59609.172414</td>\n",
       "      <td>41410.241379</td>\n",
       "      <td>42444.206897</td>\n",
       "      <td>20603.241379</td>\n",
       "      <td>9849.400000</td>\n",
       "      <td>중곡2동</td>\n",
       "      <td>신성시장앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2</td>\n",
       "      <td>101000252</td>\n",
       "      <td>5118</td>\n",
       "      <td>2703.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>11004.166667</td>\n",
       "      <td>28241.60000</td>\n",
       "      <td>45374.814815</td>\n",
       "      <td>...</td>\n",
       "      <td>69060.000000</td>\n",
       "      <td>70571.666667</td>\n",
       "      <td>76512.074074</td>\n",
       "      <td>59903.518519</td>\n",
       "      <td>40178.148148</td>\n",
       "      <td>36353.333333</td>\n",
       "      <td>23193.000000</td>\n",
       "      <td>9068.076923</td>\n",
       "      <td>중곡2동</td>\n",
       "      <td>신성시장앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3</td>\n",
       "      <td>101000252</td>\n",
       "      <td>5118</td>\n",
       "      <td>3681.083333</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>11337.857143</td>\n",
       "      <td>22246.87500</td>\n",
       "      <td>45360.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>56242.636364</td>\n",
       "      <td>52605.090909</td>\n",
       "      <td>62490.954545</td>\n",
       "      <td>49755.000000</td>\n",
       "      <td>32210.666667</td>\n",
       "      <td>24723.360000</td>\n",
       "      <td>18744.360000</td>\n",
       "      <td>8492.166667</td>\n",
       "      <td>중곡2동</td>\n",
       "      <td>신성시장앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>4</td>\n",
       "      <td>101000252</td>\n",
       "      <td>5118</td>\n",
       "      <td>5429.600000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>11256.428571</td>\n",
       "      <td>23607.00000</td>\n",
       "      <td>55623.500000</td>\n",
       "      <td>...</td>\n",
       "      <td>49254.545455</td>\n",
       "      <td>49246.521739</td>\n",
       "      <td>50751.130435</td>\n",
       "      <td>44904.521739</td>\n",
       "      <td>37905.500000</td>\n",
       "      <td>28352.680000</td>\n",
       "      <td>20732.884615</td>\n",
       "      <td>10330.560000</td>\n",
       "      <td>중곡2동</td>\n",
       "      <td>신성시장앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>5</td>\n",
       "      <td>101000252</td>\n",
       "      <td>5118</td>\n",
       "      <td>2350.636364</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>5432.500000</td>\n",
       "      <td>11010.47619</td>\n",
       "      <td>25561.904762</td>\n",
       "      <td>...</td>\n",
       "      <td>26416.090909</td>\n",
       "      <td>25811.318182</td>\n",
       "      <td>26384.454545</td>\n",
       "      <td>21789.181818</td>\n",
       "      <td>21223.583333</td>\n",
       "      <td>13381.333333</td>\n",
       "      <td>9138.260870</td>\n",
       "      <td>3998.590909</td>\n",
       "      <td>중곡2동</td>\n",
       "      <td>신성시장앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3223</th>\n",
       "      <td>8</td>\n",
       "      <td>104900129</td>\n",
       "      <td>5999</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.00000</td>\n",
       "      <td>10588.846154</td>\n",
       "      <td>...</td>\n",
       "      <td>22800.125000</td>\n",
       "      <td>29445.080000</td>\n",
       "      <td>35284.760000</td>\n",
       "      <td>30714.642857</td>\n",
       "      <td>18086.400000</td>\n",
       "      <td>14330.960000</td>\n",
       "      <td>14714.346154</td>\n",
       "      <td>5100.925926</td>\n",
       "      <td>구의3동</td>\n",
       "      <td>강변역.테크노마트앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3224</th>\n",
       "      <td>9</td>\n",
       "      <td>104900129</td>\n",
       "      <td>5999</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.00000</td>\n",
       "      <td>11197.333333</td>\n",
       "      <td>...</td>\n",
       "      <td>22315.068966</td>\n",
       "      <td>30135.724138</td>\n",
       "      <td>38724.827586</td>\n",
       "      <td>32552.000000</td>\n",
       "      <td>25597.400000</td>\n",
       "      <td>18897.666667</td>\n",
       "      <td>17745.700000</td>\n",
       "      <td>6152.533333</td>\n",
       "      <td>구의3동</td>\n",
       "      <td>강변역.테크노마트앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3225</th>\n",
       "      <td>10</td>\n",
       "      <td>104900129</td>\n",
       "      <td>5999</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.00000</td>\n",
       "      <td>16298.258065</td>\n",
       "      <td>...</td>\n",
       "      <td>34483.387097</td>\n",
       "      <td>42710.419355</td>\n",
       "      <td>70004.967742</td>\n",
       "      <td>50100.000000</td>\n",
       "      <td>37004.800000</td>\n",
       "      <td>28757.419355</td>\n",
       "      <td>28204.838710</td>\n",
       "      <td>9611.032258</td>\n",
       "      <td>구의3동</td>\n",
       "      <td>강변역.테크노마트앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3226</th>\n",
       "      <td>11</td>\n",
       "      <td>104900129</td>\n",
       "      <td>5999</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.00000</td>\n",
       "      <td>30937.066667</td>\n",
       "      <td>...</td>\n",
       "      <td>58794.571429</td>\n",
       "      <td>71392.714286</td>\n",
       "      <td>108013.714286</td>\n",
       "      <td>82263.310345</td>\n",
       "      <td>55120.689655</td>\n",
       "      <td>45304.413793</td>\n",
       "      <td>39824.931034</td>\n",
       "      <td>14025.103448</td>\n",
       "      <td>구의3동</td>\n",
       "      <td>강변역.테크노마트앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3227</th>\n",
       "      <td>12</td>\n",
       "      <td>104900129</td>\n",
       "      <td>5999</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.00000</td>\n",
       "      <td>25904.000000</td>\n",
       "      <td>...</td>\n",
       "      <td>51528.888889</td>\n",
       "      <td>69076.800000</td>\n",
       "      <td>95540.000000</td>\n",
       "      <td>68944.800000</td>\n",
       "      <td>50037.900000</td>\n",
       "      <td>38784.600000</td>\n",
       "      <td>32420.500000</td>\n",
       "      <td>12043.500000</td>\n",
       "      <td>구의3동</td>\n",
       "      <td>강변역.테크노마트앞</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3228 rows × 29 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      사용년월  표준버스정류장ID  버스정류장ARS번호            0    1    2    3             4  \\\n",
       "0        1  101000252        5118  3944.137931  0.0  0.0  0.0  11022.814815   \n",
       "1        2  101000252        5118  2703.000000  0.0  0.0  0.0  11004.166667   \n",
       "2        3  101000252        5118  3681.083333  0.0  0.0  0.0  11337.857143   \n",
       "3        4  101000252        5118  5429.600000  0.0  0.0  0.0  11256.428571   \n",
       "4        5  101000252        5118  2350.636364  0.0  0.0  0.0   5432.500000   \n",
       "...    ...        ...         ...          ...  ...  ...  ...           ...   \n",
       "3223     8  104900129        5999     0.000000  0.0  0.0  0.0      0.000000   \n",
       "3224     9  104900129        5999     0.000000  0.0  0.0  0.0      0.000000   \n",
       "3225    10  104900129        5999     0.000000  0.0  0.0  0.0      0.000000   \n",
       "3226    11  104900129        5999     0.000000  0.0  0.0  0.0      0.000000   \n",
       "3227    12  104900129        5999     0.000000  0.0  0.0  0.0      0.000000   \n",
       "\n",
       "                5             6  ...            16            17  \\\n",
       "0     28014.00000  49520.413793  ...  73242.620690  81659.310345   \n",
       "1     28241.60000  45374.814815  ...  69060.000000  70571.666667   \n",
       "2     22246.87500  45360.000000  ...  56242.636364  52605.090909   \n",
       "3     23607.00000  55623.500000  ...  49254.545455  49246.521739   \n",
       "4     11010.47619  25561.904762  ...  26416.090909  25811.318182   \n",
       "...           ...           ...  ...           ...           ...   \n",
       "3223      0.00000  10588.846154  ...  22800.125000  29445.080000   \n",
       "3224      0.00000  11197.333333  ...  22315.068966  30135.724138   \n",
       "3225      0.00000  16298.258065  ...  34483.387097  42710.419355   \n",
       "3226      0.00000  30937.066667  ...  58794.571429  71392.714286   \n",
       "3227      0.00000  25904.000000  ...  51528.888889  69076.800000   \n",
       "\n",
       "                 18            19            20            21            22  \\\n",
       "0      76253.965517  59609.172414  41410.241379  42444.206897  20603.241379   \n",
       "1      76512.074074  59903.518519  40178.148148  36353.333333  23193.000000   \n",
       "2      62490.954545  49755.000000  32210.666667  24723.360000  18744.360000   \n",
       "3      50751.130435  44904.521739  37905.500000  28352.680000  20732.884615   \n",
       "4      26384.454545  21789.181818  21223.583333  13381.333333   9138.260870   \n",
       "...             ...           ...           ...           ...           ...   \n",
       "3223   35284.760000  30714.642857  18086.400000  14330.960000  14714.346154   \n",
       "3224   38724.827586  32552.000000  25597.400000  18897.666667  17745.700000   \n",
       "3225   70004.967742  50100.000000  37004.800000  28757.419355  28204.838710   \n",
       "3226  108013.714286  82263.310345  55120.689655  45304.413793  39824.931034   \n",
       "3227   95540.000000  68944.800000  50037.900000  38784.600000  32420.500000   \n",
       "\n",
       "                23   행정동          역명  \n",
       "0      9849.400000  중곡2동       신성시장앞  \n",
       "1      9068.076923  중곡2동       신성시장앞  \n",
       "2      8492.166667  중곡2동       신성시장앞  \n",
       "3     10330.560000  중곡2동       신성시장앞  \n",
       "4      3998.590909  중곡2동       신성시장앞  \n",
       "...            ...   ...         ...  \n",
       "3223   5100.925926  구의3동  강변역.테크노마트앞  \n",
       "3224   6152.533333  구의3동  강변역.테크노마트앞  \n",
       "3225   9611.032258  구의3동  강변역.테크노마트앞  \n",
       "3226  14025.103448  구의3동  강변역.테크노마트앞  \n",
       "3227  12043.500000  구의3동  강변역.테크노마트앞  \n",
       "\n",
       "[3228 rows x 29 columns]"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bus_dust2_5"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "ed731b4c",
   "metadata": {},
   "outputs": [],
   "source": [
    "bus_dust2_5.to_csv('bus_dust2_5.csv')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "05f92b38",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "094873c5",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e5db5b44",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "faca8f00",
   "metadata": {},
   "source": [
    "## 승하차수와 (초)미세먼지 농도를 곱한 승하차수 데이터 비교하기"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "60395b17",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>표준버스정류장ID</th>\n",
       "      <th>역명</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>101000252</td>\n",
       "      <td>신성시장앞</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>104000001</td>\n",
       "      <td>군자교입구</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>104000002</td>\n",
       "      <td>군자교입구</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>104000003</td>\n",
       "      <td>군자역.용마초등학교</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>104000004</td>\n",
       "      <td>군자역.용마초등학교</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>264</th>\n",
       "      <td>104900125</td>\n",
       "      <td>중곡119안전센터</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>265</th>\n",
       "      <td>104900126</td>\n",
       "      <td>광장자이아파트</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>266</th>\n",
       "      <td>104900127</td>\n",
       "      <td>광장자이아파트</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>267</th>\n",
       "      <td>104900128</td>\n",
       "      <td>중곡아파트.기점</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>268</th>\n",
       "      <td>104900129</td>\n",
       "      <td>강변역.테크노마트앞</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>269 rows × 2 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "     표준버스정류장ID          역명\n",
       "0    101000252       신성시장앞\n",
       "1    104000001       군자교입구\n",
       "2    104000002       군자교입구\n",
       "3    104000003  군자역.용마초등학교\n",
       "4    104000004  군자역.용마초등학교\n",
       "..         ...         ...\n",
       "264  104900125   중곡119안전센터\n",
       "265  104900126     광장자이아파트\n",
       "266  104900127     광장자이아파트\n",
       "267  104900128    중곡아파트.기점\n",
       "268  104900129  강변역.테크노마트앞\n",
       "\n",
       "[269 rows x 2 columns]"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bus_table = pd.read_csv('busdata.csv')\n",
    "bus_table = bus_table.loc[:,['표준버스정류장ID','역명']]\n",
    "bus_table"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "id": "3d276bff",
   "metadata": {},
   "outputs": [],
   "source": [
    "bususer_GJ['시간합'] = bususer_GJ.iloc[:,3:27].sum(axis=1)\n",
    "bususer_GJ\n",
    "bususer_GJ_rank = bususer_GJ[['표준버스정류장ID', '시간합']].groupby(['표준버스정류장ID']).mean().reset_index()\n",
    "bususer_GJ_rank = pd.merge(bususer_GJ_rank, bus_table, how='outer')\n",
    "bususer_GJ_rank = pd.merge(bususer_GJ_rank, bususer_GJ.loc[:,['표준버스정류장ID', '행정동']].drop_duplicates(), how='outer')\n",
    "bususer_GJ_rank = bususer_GJ_rank.sort_values(by='시간합', ascending = False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "id": "0aa17c54",
   "metadata": {},
   "outputs": [],
   "source": [
    "bus_dust2_5['시간합'] = bus_dust2_5.iloc[:,3:27].sum(axis=1)\n",
    "bus_dust2_5\n",
    "bus_dust2_5_rank = bus_dust2_5[['표준버스정류장ID', '시간합']].groupby(['표준버스정류장ID']).mean().reset_index()\n",
    "bus_dust2_5_rank = pd.merge(bus_dust2_5_rank, bus_table, how='outer')\n",
    "bus_dust2_5_rank = pd.merge(bus_dust2_5_rank,bususer_GJ.loc[:,['표준버스정류장ID', '행정동']].drop_duplicates(), how='outer')\n",
    "bus_dust2_5_rank = bus_dust2_5_rank.sort_values(by='시간합', ascending = False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "id": "a2c20822",
   "metadata": {},
   "outputs": [],
   "source": [
    "bus_dust10['시간합'] = bus_dust10.iloc[:,3:27].sum(axis=1)\n",
    "bus_dust10\n",
    "bus_dust10_rank = bus_dust10[['표준버스정류장ID', '시간합']].groupby(['표준버스정류장ID']).mean().reset_index()\n",
    "bus_dust10_rank = pd.merge(bus_dust10_rank, bus_table, how='outer')\n",
    "bus_dust10_rank = pd.merge(bus_dust10_rank, bususer_GJ.loc[:,['표준버스정류장ID', '행정동']].drop_duplicates(), how='outer')\n",
    "bus_dust10_rank = bus_dust10_rank.sort_values(by='시간합', ascending = False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "id": "1ccdeeb9",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>표준버스정류장ID</th>\n",
       "      <th>시간합</th>\n",
       "      <th>역명</th>\n",
       "      <th>행정동</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>114</th>\n",
       "      <td>104000139</td>\n",
       "      <td>147535.166667</td>\n",
       "      <td>건대입구역사거리.건대병원</td>\n",
       "      <td>화양동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>104000007</td>\n",
       "      <td>146950.166667</td>\n",
       "      <td>어린이대공원후문아차산역</td>\n",
       "      <td>능동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>100</th>\n",
       "      <td>104000125</td>\n",
       "      <td>115582.500000</td>\n",
       "      <td>건대로데오거리입구</td>\n",
       "      <td>자양4동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>104000006</td>\n",
       "      <td>102907.333333</td>\n",
       "      <td>어린이대공원후문아차산역</td>\n",
       "      <td>구의2동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>104000004</td>\n",
       "      <td>101130.833333</td>\n",
       "      <td>군자역.용마초등학교</td>\n",
       "      <td>능동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>53</th>\n",
       "      <td>104000055</td>\n",
       "      <td>95683.666667</td>\n",
       "      <td>광진구청</td>\n",
       "      <td>구의1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>89</th>\n",
       "      <td>104000111</td>\n",
       "      <td>92268.250000</td>\n",
       "      <td>광진문화예술회관</td>\n",
       "      <td>자양3동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>69</th>\n",
       "      <td>104000082</td>\n",
       "      <td>87388.333333</td>\n",
       "      <td>광진구청</td>\n",
       "      <td>자양1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>104000001</td>\n",
       "      <td>83257.000000</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>중곡1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>104000003</td>\n",
       "      <td>79825.250000</td>\n",
       "      <td>군자역.용마초등학교</td>\n",
       "      <td>중곡2동</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     표준버스정류장ID            시간합             역명   행정동\n",
       "114  104000139  147535.166667  건대입구역사거리.건대병원   화양동\n",
       "7    104000007  146950.166667   어린이대공원후문아차산역    능동\n",
       "100  104000125  115582.500000      건대로데오거리입구  자양4동\n",
       "6    104000006  102907.333333   어린이대공원후문아차산역  구의2동\n",
       "4    104000004  101130.833333     군자역.용마초등학교    능동\n",
       "53   104000055   95683.666667           광진구청  구의1동\n",
       "89   104000111   92268.250000       광진문화예술회관  자양3동\n",
       "69   104000082   87388.333333           광진구청  자양1동\n",
       "1    104000001   83257.000000          군자교입구  중곡1동\n",
       "3    104000003   79825.250000     군자역.용마초등학교  중곡2동"
      ]
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bususer_GJ_rank.head(10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "id": "66b42bdd",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>표준버스정류장ID</th>\n",
       "      <th>시간합</th>\n",
       "      <th>역명</th>\n",
       "      <th>행정동</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>89</th>\n",
       "      <td>104000111</td>\n",
       "      <td>1.748432e+06</td>\n",
       "      <td>광진문화예술회관</td>\n",
       "      <td>자양3동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>53</th>\n",
       "      <td>104000055</td>\n",
       "      <td>1.733421e+06</td>\n",
       "      <td>광진구청</td>\n",
       "      <td>구의1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>104000007</td>\n",
       "      <td>1.675926e+06</td>\n",
       "      <td>어린이대공원후문아차산역</td>\n",
       "      <td>능동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>114</th>\n",
       "      <td>104000139</td>\n",
       "      <td>1.539026e+06</td>\n",
       "      <td>건대입구역사거리.건대병원</td>\n",
       "      <td>화양동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50</th>\n",
       "      <td>104000052</td>\n",
       "      <td>1.401641e+06</td>\n",
       "      <td>광진경찰서</td>\n",
       "      <td>구의1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>104000006</td>\n",
       "      <td>1.352889e+06</td>\n",
       "      <td>어린이대공원후문아차산역</td>\n",
       "      <td>구의2동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>104000003</td>\n",
       "      <td>1.323399e+06</td>\n",
       "      <td>군자역.용마초등학교</td>\n",
       "      <td>중곡2동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>100</th>\n",
       "      <td>104000125</td>\n",
       "      <td>1.320687e+06</td>\n",
       "      <td>건대로데오거리입구</td>\n",
       "      <td>자양4동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>104000002</td>\n",
       "      <td>1.311136e+06</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>군자동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>69</th>\n",
       "      <td>104000082</td>\n",
       "      <td>1.264530e+06</td>\n",
       "      <td>광진구청</td>\n",
       "      <td>자양1동</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     표준버스정류장ID           시간합             역명   행정동\n",
       "89   104000111  1.748432e+06       광진문화예술회관  자양3동\n",
       "53   104000055  1.733421e+06           광진구청  구의1동\n",
       "7    104000007  1.675926e+06   어린이대공원후문아차산역    능동\n",
       "114  104000139  1.539026e+06  건대입구역사거리.건대병원   화양동\n",
       "50   104000052  1.401641e+06          광진경찰서  구의1동\n",
       "6    104000006  1.352889e+06   어린이대공원후문아차산역  구의2동\n",
       "3    104000003  1.323399e+06     군자역.용마초등학교  중곡2동\n",
       "100  104000125  1.320687e+06      건대로데오거리입구  자양4동\n",
       "2    104000002  1.311136e+06          군자교입구   군자동\n",
       "69   104000082  1.264530e+06           광진구청  자양1동"
      ]
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bus_dust2_5_rank.head(10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "id": "d211d8a6",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>표준버스정류장ID</th>\n",
       "      <th>시간합</th>\n",
       "      <th>역명</th>\n",
       "      <th>행정동</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>53</th>\n",
       "      <td>104000055</td>\n",
       "      <td>2.809866e+06</td>\n",
       "      <td>광진구청</td>\n",
       "      <td>구의1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>114</th>\n",
       "      <td>104000139</td>\n",
       "      <td>2.560465e+06</td>\n",
       "      <td>건대입구역사거리.건대병원</td>\n",
       "      <td>화양동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>89</th>\n",
       "      <td>104000111</td>\n",
       "      <td>2.273131e+06</td>\n",
       "      <td>광진문화예술회관</td>\n",
       "      <td>자양3동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50</th>\n",
       "      <td>104000052</td>\n",
       "      <td>2.261989e+06</td>\n",
       "      <td>광진경찰서</td>\n",
       "      <td>구의1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>104000007</td>\n",
       "      <td>1.985874e+06</td>\n",
       "      <td>어린이대공원후문아차산역</td>\n",
       "      <td>능동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>104000006</td>\n",
       "      <td>1.742298e+06</td>\n",
       "      <td>어린이대공원후문아차산역</td>\n",
       "      <td>구의2동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>104000002</td>\n",
       "      <td>1.670410e+06</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>군자동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>69</th>\n",
       "      <td>104000082</td>\n",
       "      <td>1.627444e+06</td>\n",
       "      <td>광진구청</td>\n",
       "      <td>자양1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>104000003</td>\n",
       "      <td>1.624308e+06</td>\n",
       "      <td>군자역.용마초등학교</td>\n",
       "      <td>중곡2동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>133</th>\n",
       "      <td>104000160</td>\n",
       "      <td>1.583756e+06</td>\n",
       "      <td>능동사거리.군자역</td>\n",
       "      <td>군자동</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     표준버스정류장ID           시간합             역명   행정동\n",
       "53   104000055  2.809866e+06           광진구청  구의1동\n",
       "114  104000139  2.560465e+06  건대입구역사거리.건대병원   화양동\n",
       "89   104000111  2.273131e+06       광진문화예술회관  자양3동\n",
       "50   104000052  2.261989e+06          광진경찰서  구의1동\n",
       "7    104000007  1.985874e+06   어린이대공원후문아차산역    능동\n",
       "6    104000006  1.742298e+06   어린이대공원후문아차산역  구의2동\n",
       "2    104000002  1.670410e+06          군자교입구   군자동\n",
       "69   104000082  1.627444e+06           광진구청  자양1동\n",
       "3    104000003  1.624308e+06     군자역.용마초등학교  중곡2동\n",
       "133  104000160  1.583756e+06      능동사거리.군자역   군자동"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bus_dust10_rank.head(10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4402abc8",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "id": "6438c2a1",
   "metadata": {},
   "source": [
    "## 미세먼지x인구밀집도 와 초미세먼지x인구밀집도 단순합"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "id": "9fad01ee",
   "metadata": {},
   "outputs": [],
   "source": [
    "bus_dust = bus_dust10.loc[:,['시간합']]+bus_dust2_5.loc[:,['시간합']]\n",
    "bus_dust['표준버스정류장ID'] = bususer_GJ['표준버스정류장ID']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "id": "110c6c51",
   "metadata": {},
   "outputs": [],
   "source": [
    "bus_dust_rank = bus_dust[['표준버스정류장ID', '시간합']].groupby(['표준버스정류장ID']).mean().reset_index()\n",
    "bus_dust_rank = pd.merge(bus_dust_rank,bus_table, how='outer')\n",
    "bus_dust_rank = pd.merge(bus_dust_rank,bususer_GJ.loc[:,['표준버스정류장ID', '행정동']].drop_duplicates(), how='outer')\n",
    "bus_dust_rank = bus_dust_rank.sort_values(by='시간합', ascending = False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "id": "3849c6a9",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>표준버스정류장ID</th>\n",
       "      <th>시간합</th>\n",
       "      <th>역명</th>\n",
       "      <th>행정동</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>53</th>\n",
       "      <td>104000055</td>\n",
       "      <td>4.543287e+06</td>\n",
       "      <td>광진구청</td>\n",
       "      <td>구의1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>114</th>\n",
       "      <td>104000139</td>\n",
       "      <td>4.099491e+06</td>\n",
       "      <td>건대입구역사거리.건대병원</td>\n",
       "      <td>화양동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>89</th>\n",
       "      <td>104000111</td>\n",
       "      <td>4.021563e+06</td>\n",
       "      <td>광진문화예술회관</td>\n",
       "      <td>자양3동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50</th>\n",
       "      <td>104000052</td>\n",
       "      <td>3.663630e+06</td>\n",
       "      <td>광진경찰서</td>\n",
       "      <td>구의1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>104000007</td>\n",
       "      <td>3.661800e+06</td>\n",
       "      <td>어린이대공원후문아차산역</td>\n",
       "      <td>능동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>104000006</td>\n",
       "      <td>3.095187e+06</td>\n",
       "      <td>어린이대공원후문아차산역</td>\n",
       "      <td>구의2동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>104000002</td>\n",
       "      <td>2.981546e+06</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>군자동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>104000003</td>\n",
       "      <td>2.947707e+06</td>\n",
       "      <td>군자역.용마초등학교</td>\n",
       "      <td>중곡2동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>100</th>\n",
       "      <td>104000125</td>\n",
       "      <td>2.896831e+06</td>\n",
       "      <td>건대로데오거리입구</td>\n",
       "      <td>자양4동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>69</th>\n",
       "      <td>104000082</td>\n",
       "      <td>2.891974e+06</td>\n",
       "      <td>광진구청</td>\n",
       "      <td>자양1동</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     표준버스정류장ID           시간합             역명   행정동\n",
       "53   104000055  4.543287e+06           광진구청  구의1동\n",
       "114  104000139  4.099491e+06  건대입구역사거리.건대병원   화양동\n",
       "89   104000111  4.021563e+06       광진문화예술회관  자양3동\n",
       "50   104000052  3.663630e+06          광진경찰서  구의1동\n",
       "7    104000007  3.661800e+06   어린이대공원후문아차산역    능동\n",
       "6    104000006  3.095187e+06   어린이대공원후문아차산역  구의2동\n",
       "2    104000002  2.981546e+06          군자교입구   군자동\n",
       "3    104000003  2.947707e+06     군자역.용마초등학교  중곡2동\n",
       "100  104000125  2.896831e+06      건대로데오거리입구  자양4동\n",
       "69   104000082  2.891974e+06           광진구청  자양1동"
      ]
     },
     "execution_count": 30,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bus_dust_rank.head(10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "id": "9dbaea75",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>표준버스정류장ID</th>\n",
       "      <th>시간합</th>\n",
       "      <th>역명</th>\n",
       "      <th>행정동</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>114</th>\n",
       "      <td>104000139</td>\n",
       "      <td>147535.166667</td>\n",
       "      <td>건대입구역사거리.건대병원</td>\n",
       "      <td>화양동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>104000007</td>\n",
       "      <td>146950.166667</td>\n",
       "      <td>어린이대공원후문아차산역</td>\n",
       "      <td>능동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>100</th>\n",
       "      <td>104000125</td>\n",
       "      <td>115582.500000</td>\n",
       "      <td>건대로데오거리입구</td>\n",
       "      <td>자양4동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>104000006</td>\n",
       "      <td>102907.333333</td>\n",
       "      <td>어린이대공원후문아차산역</td>\n",
       "      <td>구의2동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>104000004</td>\n",
       "      <td>101130.833333</td>\n",
       "      <td>군자역.용마초등학교</td>\n",
       "      <td>능동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>53</th>\n",
       "      <td>104000055</td>\n",
       "      <td>95683.666667</td>\n",
       "      <td>광진구청</td>\n",
       "      <td>구의1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>89</th>\n",
       "      <td>104000111</td>\n",
       "      <td>92268.250000</td>\n",
       "      <td>광진문화예술회관</td>\n",
       "      <td>자양3동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>69</th>\n",
       "      <td>104000082</td>\n",
       "      <td>87388.333333</td>\n",
       "      <td>광진구청</td>\n",
       "      <td>자양1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>104000001</td>\n",
       "      <td>83257.000000</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>중곡1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>104000003</td>\n",
       "      <td>79825.250000</td>\n",
       "      <td>군자역.용마초등학교</td>\n",
       "      <td>중곡2동</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     표준버스정류장ID            시간합             역명   행정동\n",
       "114  104000139  147535.166667  건대입구역사거리.건대병원   화양동\n",
       "7    104000007  146950.166667   어린이대공원후문아차산역    능동\n",
       "100  104000125  115582.500000      건대로데오거리입구  자양4동\n",
       "6    104000006  102907.333333   어린이대공원후문아차산역  구의2동\n",
       "4    104000004  101130.833333     군자역.용마초등학교    능동\n",
       "53   104000055   95683.666667           광진구청  구의1동\n",
       "89   104000111   92268.250000       광진문화예술회관  자양3동\n",
       "69   104000082   87388.333333           광진구청  자양1동\n",
       "1    104000001   83257.000000          군자교입구  중곡1동\n",
       "3    104000003   79825.250000     군자역.용마초등학교  중곡2동"
      ]
     },
     "execution_count": 31,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bususer_GJ_rank.head(10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "id": "9309afc7",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>표준버스정류장ID</th>\n",
       "      <th>시간합</th>\n",
       "      <th>역명</th>\n",
       "      <th>행정동</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>53</th>\n",
       "      <td>104000055</td>\n",
       "      <td>4.543287e+06</td>\n",
       "      <td>광진구청</td>\n",
       "      <td>구의1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>114</th>\n",
       "      <td>104000139</td>\n",
       "      <td>4.099491e+06</td>\n",
       "      <td>건대입구역사거리.건대병원</td>\n",
       "      <td>화양동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>89</th>\n",
       "      <td>104000111</td>\n",
       "      <td>4.021563e+06</td>\n",
       "      <td>광진문화예술회관</td>\n",
       "      <td>자양3동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50</th>\n",
       "      <td>104000052</td>\n",
       "      <td>3.663630e+06</td>\n",
       "      <td>광진경찰서</td>\n",
       "      <td>구의1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>104000007</td>\n",
       "      <td>3.661800e+06</td>\n",
       "      <td>어린이대공원후문아차산역</td>\n",
       "      <td>능동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>104000006</td>\n",
       "      <td>3.095187e+06</td>\n",
       "      <td>어린이대공원후문아차산역</td>\n",
       "      <td>구의2동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>104000002</td>\n",
       "      <td>2.981546e+06</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>군자동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>104000003</td>\n",
       "      <td>2.947707e+06</td>\n",
       "      <td>군자역.용마초등학교</td>\n",
       "      <td>중곡2동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>100</th>\n",
       "      <td>104000125</td>\n",
       "      <td>2.896831e+06</td>\n",
       "      <td>건대로데오거리입구</td>\n",
       "      <td>자양4동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>69</th>\n",
       "      <td>104000082</td>\n",
       "      <td>2.891974e+06</td>\n",
       "      <td>광진구청</td>\n",
       "      <td>자양1동</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     표준버스정류장ID           시간합             역명   행정동\n",
       "53   104000055  4.543287e+06           광진구청  구의1동\n",
       "114  104000139  4.099491e+06  건대입구역사거리.건대병원   화양동\n",
       "89   104000111  4.021563e+06       광진문화예술회관  자양3동\n",
       "50   104000052  3.663630e+06          광진경찰서  구의1동\n",
       "7    104000007  3.661800e+06   어린이대공원후문아차산역    능동\n",
       "6    104000006  3.095187e+06   어린이대공원후문아차산역  구의2동\n",
       "2    104000002  2.981546e+06          군자교입구   군자동\n",
       "3    104000003  2.947707e+06     군자역.용마초등학교  중곡2동\n",
       "100  104000125  2.896831e+06      건대로데오거리입구  자양4동\n",
       "69   104000082  2.891974e+06           광진구청  자양1동"
      ]
     },
     "execution_count": 32,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bus_dust_rank.head(10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "id": "28adc70a",
   "metadata": {},
   "outputs": [],
   "source": [
    "bususer_GJ_rank = pd.merge(bususer_GJ_rank, bususer_GJ.loc[:,['표준버스정류장ID', '행정동']].drop_duplicates(), how='outer')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "id": "835b3458",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>표준버스정류장ID</th>\n",
       "      <th>시간합</th>\n",
       "      <th>역명</th>\n",
       "      <th>행정동</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>104000139</td>\n",
       "      <td>147535.166667</td>\n",
       "      <td>건대입구역사거리.건대병원</td>\n",
       "      <td>화양동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>104000007</td>\n",
       "      <td>146950.166667</td>\n",
       "      <td>어린이대공원후문아차산역</td>\n",
       "      <td>능동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>104000125</td>\n",
       "      <td>115582.500000</td>\n",
       "      <td>건대로데오거리입구</td>\n",
       "      <td>자양4동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>104000006</td>\n",
       "      <td>102907.333333</td>\n",
       "      <td>어린이대공원후문아차산역</td>\n",
       "      <td>구의2동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>104000004</td>\n",
       "      <td>101130.833333</td>\n",
       "      <td>군자역.용마초등학교</td>\n",
       "      <td>능동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>264</th>\n",
       "      <td>104900057</td>\n",
       "      <td>800.166667</td>\n",
       "      <td>한강호텔입구</td>\n",
       "      <td>광장동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>265</th>\n",
       "      <td>104000275</td>\n",
       "      <td>64.500000</td>\n",
       "      <td>중곡동성당</td>\n",
       "      <td>중곡1동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>266</th>\n",
       "      <td>104000274</td>\n",
       "      <td>31.666667</td>\n",
       "      <td>종합의료복합단지</td>\n",
       "      <td>중곡3동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>267</th>\n",
       "      <td>104000272</td>\n",
       "      <td>29.333333</td>\n",
       "      <td>용곡삼거리</td>\n",
       "      <td>중곡3동</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>268</th>\n",
       "      <td>104000273</td>\n",
       "      <td>25.916667</td>\n",
       "      <td>중곡역.국립정신건강센터</td>\n",
       "      <td>중곡3동</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>269 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "     표준버스정류장ID            시간합             역명   행정동\n",
       "0    104000139  147535.166667  건대입구역사거리.건대병원   화양동\n",
       "1    104000007  146950.166667   어린이대공원후문아차산역    능동\n",
       "2    104000125  115582.500000      건대로데오거리입구  자양4동\n",
       "3    104000006  102907.333333   어린이대공원후문아차산역  구의2동\n",
       "4    104000004  101130.833333     군자역.용마초등학교    능동\n",
       "..         ...            ...            ...   ...\n",
       "264  104900057     800.166667         한강호텔입구   광장동\n",
       "265  104000275      64.500000          중곡동성당  중곡1동\n",
       "266  104000274      31.666667       종합의료복합단지  중곡3동\n",
       "267  104000272      29.333333          용곡삼거리  중곡3동\n",
       "268  104000273      25.916667   중곡역.국립정신건강센터  중곡3동\n",
       "\n",
       "[269 rows x 4 columns]"
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bususer_GJ_rank"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "id": "2324d2ad",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>사용년월</th>\n",
       "      <th>표준버스정류장ID</th>\n",
       "      <th>버스정류장ARS번호</th>\n",
       "      <th>0</th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "      <th>5</th>\n",
       "      <th>6</th>\n",
       "      <th>...</th>\n",
       "      <th>17</th>\n",
       "      <th>18</th>\n",
       "      <th>19</th>\n",
       "      <th>20</th>\n",
       "      <th>21</th>\n",
       "      <th>22</th>\n",
       "      <th>23</th>\n",
       "      <th>행정동</th>\n",
       "      <th>역명</th>\n",
       "      <th>시간합</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>1</td>\n",
       "      <td>104000002</td>\n",
       "      <td>5002</td>\n",
       "      <td>704</td>\n",
       "      <td>157</td>\n",
       "      <td>73</td>\n",
       "      <td>35</td>\n",
       "      <td>115</td>\n",
       "      <td>2043</td>\n",
       "      <td>2044</td>\n",
       "      <td>...</td>\n",
       "      <td>5664</td>\n",
       "      <td>7224</td>\n",
       "      <td>4794</td>\n",
       "      <td>3165</td>\n",
       "      <td>4063</td>\n",
       "      <td>2341</td>\n",
       "      <td>1103</td>\n",
       "      <td>군자동</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>70740</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>2</td>\n",
       "      <td>104000002</td>\n",
       "      <td>5002</td>\n",
       "      <td>565</td>\n",
       "      <td>136</td>\n",
       "      <td>70</td>\n",
       "      <td>29</td>\n",
       "      <td>97</td>\n",
       "      <td>1916</td>\n",
       "      <td>1800</td>\n",
       "      <td>...</td>\n",
       "      <td>4864</td>\n",
       "      <td>5932</td>\n",
       "      <td>4160</td>\n",
       "      <td>2761</td>\n",
       "      <td>3093</td>\n",
       "      <td>2293</td>\n",
       "      <td>1144</td>\n",
       "      <td>군자동</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>61197</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>3</td>\n",
       "      <td>104000002</td>\n",
       "      <td>5002</td>\n",
       "      <td>910</td>\n",
       "      <td>202</td>\n",
       "      <td>116</td>\n",
       "      <td>47</td>\n",
       "      <td>179</td>\n",
       "      <td>2567</td>\n",
       "      <td>1965</td>\n",
       "      <td>...</td>\n",
       "      <td>5631</td>\n",
       "      <td>6847</td>\n",
       "      <td>4672</td>\n",
       "      <td>2935</td>\n",
       "      <td>2693</td>\n",
       "      <td>2758</td>\n",
       "      <td>1914</td>\n",
       "      <td>군자동</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>70553</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>4</td>\n",
       "      <td>104000002</td>\n",
       "      <td>5002</td>\n",
       "      <td>1469</td>\n",
       "      <td>348</td>\n",
       "      <td>201</td>\n",
       "      <td>104</td>\n",
       "      <td>203</td>\n",
       "      <td>2709</td>\n",
       "      <td>2096</td>\n",
       "      <td>...</td>\n",
       "      <td>5990</td>\n",
       "      <td>7267</td>\n",
       "      <td>5356</td>\n",
       "      <td>3402</td>\n",
       "      <td>2993</td>\n",
       "      <td>3040</td>\n",
       "      <td>1962</td>\n",
       "      <td>군자동</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>77085</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>5</td>\n",
       "      <td>104000002</td>\n",
       "      <td>5002</td>\n",
       "      <td>1486</td>\n",
       "      <td>532</td>\n",
       "      <td>241</td>\n",
       "      <td>147</td>\n",
       "      <td>237</td>\n",
       "      <td>2777</td>\n",
       "      <td>2251</td>\n",
       "      <td>...</td>\n",
       "      <td>6095</td>\n",
       "      <td>7706</td>\n",
       "      <td>5583</td>\n",
       "      <td>3734</td>\n",
       "      <td>3404</td>\n",
       "      <td>3274</td>\n",
       "      <td>2325</td>\n",
       "      <td>군자동</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>82026</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>6</td>\n",
       "      <td>104000002</td>\n",
       "      <td>5002</td>\n",
       "      <td>1106</td>\n",
       "      <td>517</td>\n",
       "      <td>218</td>\n",
       "      <td>82</td>\n",
       "      <td>204</td>\n",
       "      <td>2438</td>\n",
       "      <td>2114</td>\n",
       "      <td>...</td>\n",
       "      <td>6103</td>\n",
       "      <td>7486</td>\n",
       "      <td>5565</td>\n",
       "      <td>3563</td>\n",
       "      <td>3300</td>\n",
       "      <td>3081</td>\n",
       "      <td>2157</td>\n",
       "      <td>군자동</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>78177</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>7</td>\n",
       "      <td>104000002</td>\n",
       "      <td>5002</td>\n",
       "      <td>1029</td>\n",
       "      <td>290</td>\n",
       "      <td>220</td>\n",
       "      <td>71</td>\n",
       "      <td>187</td>\n",
       "      <td>2648</td>\n",
       "      <td>2168</td>\n",
       "      <td>...</td>\n",
       "      <td>5891</td>\n",
       "      <td>7852</td>\n",
       "      <td>5529</td>\n",
       "      <td>3624</td>\n",
       "      <td>3474</td>\n",
       "      <td>3062</td>\n",
       "      <td>2005</td>\n",
       "      <td>군자동</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>79318</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>8</td>\n",
       "      <td>104000002</td>\n",
       "      <td>5002</td>\n",
       "      <td>948</td>\n",
       "      <td>206</td>\n",
       "      <td>154</td>\n",
       "      <td>62</td>\n",
       "      <td>197</td>\n",
       "      <td>2660</td>\n",
       "      <td>2131</td>\n",
       "      <td>...</td>\n",
       "      <td>5799</td>\n",
       "      <td>7459</td>\n",
       "      <td>5548</td>\n",
       "      <td>3352</td>\n",
       "      <td>3183</td>\n",
       "      <td>3009</td>\n",
       "      <td>1948</td>\n",
       "      <td>군자동</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>76881</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>32</th>\n",
       "      <td>9</td>\n",
       "      <td>104000002</td>\n",
       "      <td>5002</td>\n",
       "      <td>958</td>\n",
       "      <td>259</td>\n",
       "      <td>183</td>\n",
       "      <td>99</td>\n",
       "      <td>147</td>\n",
       "      <td>2445</td>\n",
       "      <td>2095</td>\n",
       "      <td>...</td>\n",
       "      <td>5948</td>\n",
       "      <td>7314</td>\n",
       "      <td>5604</td>\n",
       "      <td>3609</td>\n",
       "      <td>3224</td>\n",
       "      <td>3145</td>\n",
       "      <td>1946</td>\n",
       "      <td>군자동</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>78694</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33</th>\n",
       "      <td>10</td>\n",
       "      <td>104000002</td>\n",
       "      <td>5002</td>\n",
       "      <td>1101</td>\n",
       "      <td>314</td>\n",
       "      <td>237</td>\n",
       "      <td>107</td>\n",
       "      <td>171</td>\n",
       "      <td>2636</td>\n",
       "      <td>2303</td>\n",
       "      <td>...</td>\n",
       "      <td>6215</td>\n",
       "      <td>7441</td>\n",
       "      <td>5344</td>\n",
       "      <td>3688</td>\n",
       "      <td>3447</td>\n",
       "      <td>3168</td>\n",
       "      <td>2151</td>\n",
       "      <td>군자동</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>82287</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>34</th>\n",
       "      <td>11</td>\n",
       "      <td>104000002</td>\n",
       "      <td>5002</td>\n",
       "      <td>1082</td>\n",
       "      <td>255</td>\n",
       "      <td>172</td>\n",
       "      <td>72</td>\n",
       "      <td>155</td>\n",
       "      <td>2601</td>\n",
       "      <td>2202</td>\n",
       "      <td>...</td>\n",
       "      <td>6470</td>\n",
       "      <td>7549</td>\n",
       "      <td>5607</td>\n",
       "      <td>3542</td>\n",
       "      <td>3364</td>\n",
       "      <td>3145</td>\n",
       "      <td>1983</td>\n",
       "      <td>군자동</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>81387</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>35</th>\n",
       "      <td>12</td>\n",
       "      <td>104000002</td>\n",
       "      <td>5002</td>\n",
       "      <td>1010</td>\n",
       "      <td>300</td>\n",
       "      <td>171</td>\n",
       "      <td>66</td>\n",
       "      <td>160</td>\n",
       "      <td>2360</td>\n",
       "      <td>2107</td>\n",
       "      <td>...</td>\n",
       "      <td>6653</td>\n",
       "      <td>7653</td>\n",
       "      <td>5547</td>\n",
       "      <td>3515</td>\n",
       "      <td>3356</td>\n",
       "      <td>3377</td>\n",
       "      <td>2109</td>\n",
       "      <td>군자동</td>\n",
       "      <td>군자교입구</td>\n",
       "      <td>82514</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>12 rows × 30 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "    사용년월  표준버스정류장ID  버스정류장ARS번호     0    1    2    3    4     5     6  ...  \\\n",
       "24     1  104000002        5002   704  157   73   35  115  2043  2044  ...   \n",
       "25     2  104000002        5002   565  136   70   29   97  1916  1800  ...   \n",
       "26     3  104000002        5002   910  202  116   47  179  2567  1965  ...   \n",
       "27     4  104000002        5002  1469  348  201  104  203  2709  2096  ...   \n",
       "28     5  104000002        5002  1486  532  241  147  237  2777  2251  ...   \n",
       "29     6  104000002        5002  1106  517  218   82  204  2438  2114  ...   \n",
       "30     7  104000002        5002  1029  290  220   71  187  2648  2168  ...   \n",
       "31     8  104000002        5002   948  206  154   62  197  2660  2131  ...   \n",
       "32     9  104000002        5002   958  259  183   99  147  2445  2095  ...   \n",
       "33    10  104000002        5002  1101  314  237  107  171  2636  2303  ...   \n",
       "34    11  104000002        5002  1082  255  172   72  155  2601  2202  ...   \n",
       "35    12  104000002        5002  1010  300  171   66  160  2360  2107  ...   \n",
       "\n",
       "      17    18    19    20    21    22    23  행정동     역명    시간합  \n",
       "24  5664  7224  4794  3165  4063  2341  1103  군자동  군자교입구  70740  \n",
       "25  4864  5932  4160  2761  3093  2293  1144  군자동  군자교입구  61197  \n",
       "26  5631  6847  4672  2935  2693  2758  1914  군자동  군자교입구  70553  \n",
       "27  5990  7267  5356  3402  2993  3040  1962  군자동  군자교입구  77085  \n",
       "28  6095  7706  5583  3734  3404  3274  2325  군자동  군자교입구  82026  \n",
       "29  6103  7486  5565  3563  3300  3081  2157  군자동  군자교입구  78177  \n",
       "30  5891  7852  5529  3624  3474  3062  2005  군자동  군자교입구  79318  \n",
       "31  5799  7459  5548  3352  3183  3009  1948  군자동  군자교입구  76881  \n",
       "32  5948  7314  5604  3609  3224  3145  1946  군자동  군자교입구  78694  \n",
       "33  6215  7441  5344  3688  3447  3168  2151  군자동  군자교입구  82287  \n",
       "34  6470  7549  5607  3542  3364  3145  1983  군자동  군자교입구  81387  \n",
       "35  6653  7653  5547  3515  3356  3377  2109  군자동  군자교입구  82514  \n",
       "\n",
       "[12 rows x 30 columns]"
      ]
     },
     "execution_count": 35,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "bususer_GJ[bususer_GJ['표준버스정류장ID'] == 104000002]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "ba3d9afa",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4bc998f1",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "16017a75",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e0026376",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "cab430b5",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d1f262f5",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "id": "b9c59174",
   "metadata": {},
   "outputs": [],
   "source": [
    "bus_dust_top50 = bus_dust_rank.head(50)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "id": "228817e5",
   "metadata": {},
   "outputs": [],
   "source": [
    "bus_dust_top50.to_csv('bus_dust_top50.csv')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "33504eb6",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "aa947a43",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "id": "9f6eaa66",
   "metadata": {},
   "outputs": [],
   "source": [
    "new4_2 = bus_dust10.copy()\n",
    "\n",
    "new4_2.loc['시간 별 평균 이용자 수',:]= new4_2.loc[:,'0':'23'].mean(axis=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "id": "268eb5e3",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x24b77ed2c70>]"
      ]
     },
     "execution_count": 59,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3dd3xc1Znw8d+jbnXJKh5LsuUqW7blItkQIHRigk0LzWx2w24KuwQCIRsSeNklyZtl32Rhk8BmQ8ImJCFLAAcwppqWUJKAsVzV3JuqJVlWsaw6c94/5o4ZhOr00Tzfz2c+ujpz59xn7Ktn7pxz7jlijEEppVRkiAp2AEoppQJHk75SSkUQTfpKKRVBNOkrpVQE0aSvlFIRJCbYAYwlKyvLFBYWBjsMpZQKK1u3bm01xmQPLQ/5pF9YWEh5eXmww1BKqbAiIkeGK9fmHaWUiiCa9JVSKoJo0ldKqQiiSV8ppSKIJn2llIogmvSVUiqCaNJXSqkIokk/wMoPt1FR1xHsMJRSEWrMpC8ij4lIs4hUupUtE5EPRGSHiJSLyCq35+4Rkf0iskdEVruVl4pIhfXcwyIivn87oe9bz+zi3ucrgh2GUipCjedK/zfApUPK/gP4njFmGXCf9TsiUgysAxZZr/mZiERbr3kEuBmYZz2G1jnpdfcNcuh4N1UNnXT3DQY7HKVUBBoz6Rtj3gXahhYDqdZ2GtBgbV8JPGWM6TPGHAL2A6tExAakGmPeN86luh4HrvLFGwgnu5u6MAbsDsPO2vZgh6OUikCetul/HXhARGqBB4F7rPI8oNZtvzqrLM/aHlo+LBG52Wo2Km9pafEwxNBT3dh5envL4RNBjEQpFak8Tfq3AHcaYwqAO4FfWeXDtdObUcqHZYx51BhTZowpy87+xCRxYau6oZO0KbEsmJZC+ZGhX56UUsr/PE36NwHPWdt/AFwduXVAgdt++Tibfuqs7aHlEaWmsZNiWyplhRlsO3KCQbsj2CEppSKMp0m/ATjP2r4Q2GdtvwCsE5F4EZmFs8P2Q2NMI9AlImdao3a+AGz0Iu6wY3cYdjd1Ujw9lZWFmXT329nd1BXssJRSEWbM+fRF5EngfCBLROqA7wBfAR4SkRigF+eoHIwxVSKyHqgGBoFbjTF2q6pbcI4EmgK8aj0ixqHWbnoHHCy0pVI6MwOArUdOsDgvLciRKaUiyZhJ3xhz4whPlY6w//3A/cOUlwOLJxTdJOLqxC22pZKXPgVbWgJbDrdx01mFwQ1MKRVR9I7cAKlp7CQ2Wpibk4yIUDozg/LDJ3COYFVKqcDQpB8g1Q2dzMtJIS7G+U++sjCTps5e6tt7ghyZUiqSaNIPkOrGThbaUk//XlbobNcv1/H6SqkA0qQfAC1dfbR09VE8/aOkv2BaKsnxMWw5rOP1lVKBo0k/AGrcOnFdoqOE5TPS2XpEr/SVUoGjST8AqodJ+uBs199zrIuOnoFghKWUikCa9AOguqGTvPQppCXGfqy8bGYGxsC2o3q1r5QKDE36ATC0E9dl2Yx0oqOEcm3XV0oFiCZ9P+sdsHOw5eTHOnFdEuNiWDw9VWfcVJPe7zcf5e9+tVnvSwkBmvT9bE9TFw7zyfZ8l9KZmeysbad/UCdfU5OTMYZfvHuA9/a1skPXkQg6Tfp+5urEXTTMlT7AysIM+gYdVDbourlqctp2tJ0jx08B8PKuxiBHozTp+1l1Qycp8THkZ0wZ9vnS0zdpabu+mpw2bK8jITaKM2dn8nJFIw6HNvEEkyZ9P3N14o60DnxOSgIzpybqnblqUuofdPDSrkY+UzyNG1YW0NjRy/ZaPdeDSZO+Hzkcht2NncN24rorm5lJ+RGdfE1NPn/a00z7qQGuXpHHxQtziYuJ4sWd2sQTTJr0/eho2ym6++0jduK6rCzMoK27n4Ot3QGKTKnA2LCtnqzkOD49N4uUhFjOn5/NK9rEE1Sa9P3I1Yk73Bh9d67J17ZqE4+aRDpODfDH3c1cvnQ6MdHOVLOmxEZzVx/lOv1I0GjS96Pqhk6io4R5ucmj7jcnO5mMxFidfE1NKi9XNNJvd/C55R8tj33xwlziY6J4aVfELZEdMjTp+1FNYydzs5NJiI0edT/noiqZOvmamlQ2bK9jbk4yi/M++qabFB/DhQtyeKWiCbs28QSFJn0/qh5HJ65LWWEGB1u7aT3Z5+eolPK/2rZTbDl8gquX531i5Nrakum0nuxj86HjQYousmnS95O27n4aO3pZaEsZ1/4rdVEVNYls2F4PwFXL8z7x3AULspkSG603agWJJn0/+WgO/bRx7b84L424mCi2HtF2fRXejDFs2F7PGbMyyUv/5E2JiXExXLQwh02VTQzadfqRQNOk7yc1p0fujO9KPz4mmqX5aTr5mgp7O2rbOdTazedWfPIq32VtiY3j3f18cFAvcgJtzKQvIo+JSLOIVA4p/5qI7BGRKhH5D7fye0Rkv/XcarfyUhGpsJ57WEa6RXWSqG7oZFpqAlOT48f9mrLCTCrrO+jpt/sxMqX86/nt9cTHRPHZJbYR9zm/KIekuGgdxRME47nS/w1wqXuBiFwAXAmUGGMWAQ9a5cXAOmCR9ZqfiYhr6MojwM3APOvxsTonG+f0C+O7yncpm5nBoMOws05nIlThacDu4MVdjVxcnEtqQuyI+yXERnNxcS6bqpoY0CaegBoz6Rtj3gWGfge7BfiBMabP2qfZKr8SeMoY02eMOQTsB1aJiA1INca8b5xzDTwOXOWrNxFq+gbt7G8efg790ZTO1MnXVHh7Z08Lbd39fG6YDtyh1iyx0X5qgL/sbw1AZMrF0zb9+cCnRWSziLwjIiut8jyg1m2/Oqssz9oeWj4sEblZRMpFpLylpcXDEINn37GTDDrMuDtxXdIT45ifm6x3K6qwtWF7PZlJcZw7P3vMfc+dn01KfIyO4gkwT5N+DJABnAncBay32uiHa6c3o5QPyxjzqDGmzBhTlp099skTak4vhD7BK33g9E1aeuOKCjcdPQO8UXOMy0tsxEaPnVoSYqO5pDiX16qadBGhAPI06dcBzxmnDwEHkGWVF7jtlw80WOX5w5RPStUNnSTGRTMzM3HCr11ZmEFX7yB7j3X5ITKl/OfVikb6Bx1cvSJ/7J0ta5fa6Owd5M/7w+8bfbjyNOk/D1wIICLzgTigFXgBWCci8SIyC2eH7YfGmEagS0TOtL4RfAHY6HX0Iaq6sZMF01KIipr4AKWVhZkA2sSjws5z2+uZnZXE0vzxN2ueMzeb1IQYXtImnoAZz5DNJ4H3gSIRqRORLwGPAbOtYZxPATdZV/1VwHqgGtgE3GqMcY0/vAX4Jc7O3QPAqz5/NyHAGEPNBKZfGCo/Ywo5KfHamavCSt2JU3x4qG3YaRdGExcTxepF03ij6hi9AzpUORBixtrBGHPjCE/97Qj73w/cP0x5ObB4QtGFoboTPXT1Dk64E9dFRFhZmKnTMaiwsnGHs7V2uGkXxrKmxMYfttbx3r5WLinO9XVoagi9I9fHqid4J+5wygozqG/voaG9x1dhKeU3xhie21bHqsJMCjzoxzp7bhbpibF6o1aAaNL3seqGTqIEFkzzrHkHnMsngrbrq/BQUd/BgZZuj67yAWKjo7h00TTerNYmnkDQpO9j1Y2dzMpKYkrc6HPoj2ahLYXEuGi2aru+CgPPbasnLjqKNaNMuzCWNSU2uvvtvL2neeydlVc06fuYsxPXs/Z8l5joKFbMyNDJ11TIG7A7eHFnAxctzCEtceRpF8byqdlTyUyK01E8AaBJ34c6egaoO9Ez5kLo41E6M4PdTZ109Q74IDKl/OO9fS0c7+7nag+bdlxioqO4dPE03qpp5lT/oI+iU8PRpO9DE51OeTQrCzNxGNh+VCdfU6Frw/YGMhJjOb8ox+u61pbY6Bmw86fdeqOWP2nS96HqBs+nXxhq2Yx0oqNEx+urkNXVO8DrVU2sLZlOXIz3qeSMWVPJSo7n5QodxeNPmvR9qKaxk6zkeHJSEryuKzk+hoW2FG3XVyHr1com+gYdHo/aGSo6SrhsyTT+uLuZ7j5t4vEXTfo+NJGF0MejbGYmO2rbdb5xFZI2bKuncGoiK2ak+6zONUts9A44eGu3juLxF036PtI/6GDfsZM+ac93KSvMoGfAfrrZSKlQ0dDewweHjnPVBKddGMvKwkxyUuJ5aac28fiLJn0fOdBykn67wycjd1z0Ji0Vqp7fUY8xeD1qZ6ioKOGyJTbe3tuiI9f8RJO+j7hG7izyYfPOtLQECjKnaGeuCinGGDZsq6d0ZgYzpyb5vP7Ll9roH3TwZs0xn9etNOn7THVDJwmxUczKSvZpvWUzM9ly+ATOVSaVCr6qhk72NZ/0+VW+y/KCDGxpCbqilp9o0veR6sZOinJTiPZgDv3RlBVm0Hqyj6Ntp3xar1Ke2rC9ntho8WrahdG4mnje2dtCR4828fiaJn0fMMb4fOSOi2tRFR26qULBoN3Bxh0NXFCUQ0ZSnN+Os7bExoDd8Ea1NvH4miZ9H2js6KX91IBPO3Fd5mYnk5oQo+36KiT8eX8rrSf7+NwK/zTtuCwrSCcvfQov63TLPqdJ3wdqvFgIfSxRUUJZYaaO4FEhYcP2etKmxHLBAu+nXRiNiLC2xMZ7+1ppP9Xv12NFGk36PuAaR1/kxRz6oykrzGB/80nauvXkV8Fzsm+Q16qaWFNiIz7G86nDx2tNiY1Bh+H1Km3i8SVN+j5Q3dhJ4dREkuPHXH3SI67x+lv1al8F0WuVTfQOOPicn0btDLUkL40ZmYm8qE08PqVJ3wf81YnrUpKfRlx0FOVHtF1fBc+G7fUUZE6hdGZGQI4nIqwpsfHXA8f1W64PadL30sm+QY4cP+WXTlyXhNholuSn6WLpKmiaOnr5y4FWrl7m22kXxrK2xIbdYdhU2RSwY052mvS9tPv0HPr+S/oAZTMzqKjr0DVEVVBsdE27sCI/oMcttqUyKytJp1v2oTGTvog8JiLNIlI5zHPfFBEjIlluZfeIyH4R2SMiq93KS0WkwnruYQnk5YIfVftx5I67ssJM+u0OKuo7/HocpYazYXs9ywrSmZXl+2kXRuMaxfP+geO0dPUF9NiT1Xiu9H8DXDq0UEQKgEuAo25lxcA6YJH1mp+JiKub/xHgZmCe9fhEneGouqGTjMRYpqV6P4f+aFztqFt0vL4KsOqGTnY3dfl9bP5I1pTYcBjYVKVNPL4wZtI3xrwLDJdpfgx8C3CfFOZK4CljTJ8x5hCwH1glIjYg1RjzvnFOIvM4cJXX0YeAGqsT199fXDKT4piTncRWbddXAbZhex0xUcLakulBOX5RbgpzspN0umUf8ahNX0SuAOqNMTuHPJUH1Lr9XmeV5VnbQ8tHqv9mESkXkfKWltBdL3PQ7mB3U5dfO3HdrbRu0nI4dPI1FRh2h2HjjgbOL8om04/TLozG2cQznQ8Pt9Hc2RuUGCaTCSd9EUkE7gXuG+7pYcrMKOXDMsY8aowpM8aUZWdnTzTEgDnU2k3foMPvnbguZYWZdPQMsL/lZECOp9RfD7TS3NXH1csD24E71NoSG8Y4l2hU3vHkSn8OMAvYKSKHgXxgm4hMw3kFX+C2bz7QYJXnD1Me1gLVietSZrXr69BNFSgbttWTkhDDRQv9O+3CWOblplCUm8JLeqOW1yac9I0xFcaYHGNMoTGmEGdCX2GMaQJeANaJSLyIzMLZYfuhMaYR6BKRM61RO18ANvrubQRHdWMncdFRzMn27Rz6I5k5NZGs5HidfE0FxKn+QTZVNbFmiY2EWP9PuzCWNSU2thw+QVOHNvF4YzxDNp8E3geKRKRORL400r7GmCpgPVANbAJuNca4BpbfAvwSZ+fuAeBVL2MPuuqGTuZPSyY2OjC3O4gIZTMzdPI15Xftp/r5xTsHOdVv99tiKRO1psQ5f//LFbq4ijfGnCzGGHPjGM8XDvn9fuD+YfYrBxZPML6QZYyhuqGTC/082+BQZYUZbKpq4lhnL7l+HiaqIofDYaio7+DtPS28s7eZHbXtOIxzChDXmg7BNic7mYW2VF7e1cCXzpkV7HDCln9mCIsALV19HO/uD1h7vovrD7D88InTVz5KeaKlq4/39rXwzt4W3tvXSlt3PyJQkp/ObRfO47z52SzNTyPKx6vBeWNtiY0HXttDfXsPeelTgh1OWNKk76EqVydugEbuuBRPT2VKbDTlR9o06asJGbA72H60nXf2NvPO3hYq653ncFZyHOfPz+a8omzOmZvF1OT4IEc6MlfSf2VXI185d3awwwlLmvQ95Fo4ZWGAr/Rjo6NYVpCuI3jUuDS09/DO3hbe2dPCX/a30tU3SHSUUDojg7tWF3He/GyKbakhdTU/mplTk1iSl8ZLuxo06XtIk76Hqhs6yc+YQmpCbMCPvbIwg/9++wDdfYMk+WkOfxWeBu0OPjjYdvpqfu8x5z0dtrQE1pTYOG9+NmfNzSJtSuDPW1+5fKmNf39lNwdaTgZs5NxkohnDQ9WNnQFv2nEpLczE7tjPjtp2zp6bNfYLVEQwxvDVJ7bxevUx4qKjWDUrk+tKCzivKJt5OckBnRLZn65alscPXt3Nc9vquGv1gmCHE3Y06XvgVP8gh1q7uWJpcOYiWTEjnShxTr6mSV+5/O8HR3i9+hhfv3geX/n07En7LTAnNYFz52fz3LZ6vnFJEdFh0jQVKnQ+fQ/saerCmMB34rqkJMSyYFqqLp+oTtvT1MW/vVzD+UXZ3HHRvEmb8F2uLc2nsaOX9w8cD3YoYUeTvgeqA7RwymjKCjPYduQEg3ZH0GJQoaF3wM7tT24nJSGGB65dOmmacUZz8cJcUhNieGZr7dg7q4/RpO+B6oZOUhJiyM8I3jjhssJMuvvt7G7qCloMKjT84NXd7DnWxQPXLSU7JXSHW/pSQmw0ly+dzqaqJrp6B4IdTljRpO8BVyduMK+oVha6Jl/TeXgi2R93H+M3fz3MP5xdyAVFwZ0ULdCuLc2nd8DBKzotw4Ro0p8gu8Owp6kr4HfiDmVLm0Je+hQ2H9KkH6mau3r55h92sWBaCt++NPJGsSwrSGd2dhLPbK0be2d1mib9CTpyvJtT/fagtue7nF+UzZ/2NNOpX28jjsNh+Of1O+nuG+S/blweErNgBpqIcG1pPlsOn+Bwa3ewwwkbmvQnqDpI0y8M5/qyAnoHHLywQ+cYjzSP/eUQ7+1r5V/XFjMvNyXY4QTN1cvzEIHntunV/nhp0p+g6oZOYqKEebnBvxOwJD+NBdNSWF+uIxgiSWV9Bz/ctJtLinP5/Bkzgh1OUNnSpnDO3Cye3Vavy4iOkyb9Capu7GRuTjLxMcH/Oi0iXF9WwK66jtNzAanJ7VT/IHc8tZ3MpDh+eE1JRAzPHMu1pfnUt/fwwSEdsz8emvQnqCaI0y8M5+rlecRFR/H0Fr3ajwTff6mGg63d/Oj6ZUFbqDzUrF40jZT4GO3QHSdN+hPQerKPY519QR+54y4jKY5LFuXy/I56+gbtY79Aha1NlU08+eFRbj53tk6/4SYhNpq1S21sqmyiu28w2OGEPE36E1ATQp247m4oK6D91ACvVx0LdijKTxo7erj7uV0syUvjny8pCnY4IeeaFfmc6rfrmP1x0KQ/AdUNwZ9+YTjnzM0iL32KduhOUnaH4RtP76R/0MFD65YRF6N/tkOVzsygcGqiNvGMg549E1DT2Mn0tAQyQqwtNSrKOV75z/tbqTtxKtjhKB/7xbsHeP/gcb57xSJm6/zxw3KN2d98qI3aNv0bGI0m/QmobuwMuat8l+vK8gH4Q7le6UwmO2rb+dHre1lTYuO60vxghxPSrl6Rjwg8q2P2R6VJf5x6B+wcaOkOqU5cd/kZiZwzN4tnttZh1/HKk8LJPufwzNzUBP79qiU6PHMMeelTOGvOVJ7dVqdj9kcxZtIXkcdEpFlEKt3KHhCR3SKyS0Q2iEi623P3iMh+EdkjIqvdyktFpMJ67mEJszN477Eu7A4Tcp247q4vK6C+vYe/7G8NdijKB76zsYratlP8+IZlpCWG7/KGgXRtaT61bT1s0YkIRzSeK/3fAJcOKXsDWGyMKQH2AvcAiEgxsA5YZL3mZyLiuovpEeBmYJ71GFpnSDs9cidEr/QBPrMol/TEWJ7WDt2w98LOBp7dVsdtF85j1azMYIcTNlYvmkZSXLR26I5izKRvjHkXaBtS9roxxjUg9gPA1dh4JfCUMabPGHMI2A+sEhEbkGqMed8YY4DHgat89SYCobqhk6S4aAoyEoMdyojiY6K5alkeb1Qd40R3f7DDUR6qbTvFvRsqWDEjndsvnBvscMJKYlwMa0psvFLRyKl+HbM/HF+06X8ReNXazgPcLzPrrLI8a3to+bBE5GYRKReR8paWFh+E6D1XJ25UiK/HecPKAvrtDjZsrw92KMoDg3YHdz69A2PgoXXLiYnWbreJumZFPt39djZVNgU7lJDk1RklIvcCg8ATrqJhdjOjlA/LGPOoMabMGFOWnZ3tTYg+4XAYahqDP4f+eCy0pVKSn8b68lqcX6pUOPnpn/ZTfuQE91+9mILM0P1WGcpWFmYyIzNRR/GMwOOkLyI3AWuBz5uPsksdUOC2Wz7QYJXnD1MeFupO9HCybzCkO3HdXV9WwO6mLnbVdQQ7FDUB5YfbePitfXxueR5XLhvxi7AaQ1SUcM2KfP564LjetzIMj5K+iFwKfBu4whjj/q/6ArBOROJFZBbODtsPjTGNQJeInGmN2vkCsNHL2AOmutGZPEN1jP5QVyybTkJslHbohpGOngHueGoH+RmJfO/KRcEOJ+x9bkUexsCGbdrMOdR4hmw+CbwPFIlInYh8CfgpkAK8ISI7ROTnAMaYKmA9UA1sAm41xrhmAbsF+CXOzt0DfNQPEPKqGzqJEiiaFh6LVaQmxHLZYhsv7migp18nYQt1xhj+5flKmjp7eWjdMlISdHimtwoyEzlzdibPbqvTZs4hxjN650ZjjM0YE2uMyTfG/MoYM9cYU2CMWWY9/slt//uNMXOMMUXGmFfdysuNMYut524zYfQ/Ud3YyZzs5LBaku76lQV09Q3qBFQhzhjDb/56mBd3NvCNS+azfEZGsEOaNK5Zkc/h46fYeuREsEMJKTo0YByqGjrDohPX3RmzMimcmqhNPCGss3eA25/awfderOaComz+6bw5wQ5pUrlsiY3EuGjt0B1Ck/4YWk/20djRy5K8tGCHMiEiwnVlBXx4qI1Dumh0yNl65ASXPfQer1Q0ctfqIn5500qiQ3w4cLhJio/hs4ttvLSzUZs53WjSH0NFvbMTd3GYJX1w3pIeJeiUyyHE7jD89I/7uP4X7yMCf/inT3HrBXM14fvJNaV5dPUN8nq1jtl30aQ/hkpr2GO4Ne8A5KYmcEFRDs9urWPQ7gh2OBGvsaOHz//yAx58fS9rlth4+fZPs0Lb8P3qzFlTyUufotMyuNGkP4aK+g5mZSWRGqYjKq5fWUBzVx9v7wmNO5sj1WtVTXz2offYVdfBg9ct5aF1y8L2nAonzjH7efx5fyuNHT3BDickaNIfQ2V9R1g27bhcuCCHrOR47dANkt4BO/duqOAff7eVgoxEXvraOVxbmq/TJAfQNaX5zjH7OjUJoEl/VMdP9tHQ0cuSvPBr2nGJjY7imhV5/HF3M81dvcEOJ6Lsaeriip/+mSc2Oxczf/aWs3TlqyCYOTWJVYWZPLNVx+yDJv1RVVpr4obzlT7AdWUF2B2G5/TuxIAwxvD4+4e5/Kd/pq17gMe/uIr/c9lCXds2iK4pzeNgSzfba9uDHUrQ6Vk4isowHrnjbm5OMmUzM1i/RSdh87e27n6+8vhW7ttYxdlzprLp65/m3PnBnzQw0l22xEZCbBTPaocuMcEOIJRV1HVQODVxUnS4Xb+ygG89s4vyIydYWRh5i3I0tPfwkzf3cqi1m8V5aSzNT6ckP43CqUk+my77rwdaufPpHZzoHuC+tcX8w9mF2nYfIlISYrl00TRe3NnAv64tDqu7631Nk/4oKuo7WD4jfewdw8CaJTa+90IVT2+pjaik39U7wM/fOcAv3zuEARZNT+WpD2v59V8OA5CaEEOJ9QFQkp/OsoJ0pqUlTOgYA3YHP3lzLz97+wCzspJ47O9Xsmh6eH87nIyuLS3g+R0NvFlzjLUl04MdTtBo0h/Bie5+6tt7+LtPzQx2KD6RFB/D5Uuns3FHA9+5vHjST+o1YHfw1JZafvLGXo5393PVsul8c3UR+RmJDNod7Gs+ya66dnbWdbCrrp1H3z3IoLWYdk5KvPUBkHb6AyE9MW7Y4xw9forbn9rOjtp21q0s4L7Li0mM0z+rUPSpOVOZnpbAM1vrNOmrT3LdiRtu0y+M5vqVBTy1pZaXdjVy46oZwQ7HL4wxvFXTzP97tYYDLd2cMSuTX69ZSEn+R9/YYqKjWGhLZaEtlRtWOst6B+xUN3ayq9b5QbCzrp03a46dfs3MqYmnm4SWFqSzaHoqb1Qf494NlYjAf//NCtaU2AL9dtUEREcJV6/I45G3D3Css5fc1Il9o5ssNOmP4PT0C5Poa/rygnTm5STz9JbaSZn0K+o6uP+Vaj442Mbs7CT+5wtlXLwwZ1zt6gmx0ayYkfGxO2Q7eweorOtgR107u2o7KD/cxgs7nWv/RAk4DJTNzOAn65aRH8JrJ6uPXLMin//+0wGe317PP0boBHea9EdQWd/BjMxE0hInTzOIiHDDygL+7eUa9h7rYn5ueKwPMJb69h4efG0PG7bXk5kUx/evXMS6VTOI9XJ92dSEWM6am8VZc7NOlzV39bKr1tkklJkUx9+eOVPXsQ0js7OTWTEjnWe21nHzubMjsqNdz9YRVDZ0TKqmHZerl+cRGy08vSX879Dt7B3gh5t2c8GDb/NKRSNfPX8Ob991Pn/3qUKvE/5IclISuLg4l298poi/P3uWJvwwdG1pAfuaT57+Nh9p9IwdRvupfmrbesJ+fP5wpibHc/HCXDZsr6d/MDwnYRuwO3j8/cOc/8DbPPL2AdYusfHHb57Pty5dMCmG1yr/WlNiIz4mKmInYdOkP4zKesoRb28AABT5SURBVOeduJPxSh+cHbpt3f0f66gMB8YYXq9qYvWP3+W+jVXMz03mxdvO4Uc3LCMvfUqww1NhIm1KLJ9ZNI2NOxroG4y8efY16Q/jozn0w3fOndGcOy8bW1pCWDXx7Kpr54ZHP+Dm321FBH51UxlPfuVMluRPzg9m5V/XlubT0TPAH2uagx1KwGlH7jAq6zvIz5gy4tjscBcdJVxbms9P/7SfhvYepofwVXJ9ew8PbNrN8zsamJoUx/evWsyNKwu0LV155Zy5WeSmxvPM1jo+uySyhtrqX84wKuonZyeuu+tKCzCGkG3XNMbwzNY6Vv/4XV6tbOLWC6xOWh0to3wgOkq4enk+b+9toaWrL9jhBJT+9QzRcWqAo22nJmUnrrsZUxM5a85U1pfX4nCE1iRsbd393PK/2/jmH3ZSbEvlzW+cx12rF0z6u4hVYF1bmofdYdi4w/PZZx0OQ++Ana7eAU5099M7EPp9BGM274jIY8BaoNkYs9gqywSeBgqBw8D1xpgT1nP3AF8C7MDtxpjXrPJS4DfAFOAV4A4TglM+VjZMvjtxR3LDygLueGoH7x88ztluY9GD6U97mvnWM7toP9XP3Z9dwFc+PVvXj1V+MTcnhaUF6fzi3YN8cPA4/XbDoN3BgN1Bv90wMOhg0OFgwG7oH3SWDzqc5f3Wtn3IBdPUpDieveUsCrOSgvSuxjaeNv3fAD8FHncruxt4yxjzAxG52/r92yJSDKwDFgHTgTdFZL4xxg48AtwMfIAz6V8KvOqrN+IrlZNw+oWRrF40jdSEGJ7eUhv0pN/Tb+ffX6nhdx8cYX5uMr/9h1VhuS6xCi+3XTCX/3x9D40dvcRGRxEXHUVsdBSJcc6fcTFCTNQnt2NjhNgh21FRwsNv7eOrT2zjua+eFbIzeY6Z9I0x74pI4ZDiK4Hzre3fAm8D37bKnzLG9AGHRGQ/sEpEDgOpxpj3AUTkceAqQjDpV9R3kJc+hYykydmJ6y4hNpqrlufx1JZaOk4NBO3u45217dz59A4Otnbz5XNm8c3VRSH7B6Mml0uKc7mkONdn9c3MTOTLj5fz/Zequf/qJT6r15c8bdPPNcY0Alg/c6zyPMB9HGCdVZZnbQ8tH5aI3Cwi5SJS3tIS2AW9KyOgE9fd9WUF9A86eN6Ldk1PDdodPPTmPj73yF/pGbDz+y+fwb9E+FznKrxdXJzLzefO5onNR73qK/AnX3fkDtf4akYpH5Yx5lFjTJkxpiw7O3CrDnX2DnD4+KmIGvu9OC+NRdNTAz5m/1BrN9f+/H1+/OZe1pbY2HTHuR+b40apcHXX6iJKZ2bwf56r4EDLyWCH8wmeJv1jImIDsH667nCoAwrc9ssHGqzy/GHKQ4qrPX9RhLUl37CygOrGztPv35+MMTyx+QiXPfQeB1tO8vCNy3lo3fJJNbGdimyx0VH8143LiYuJ4tYntoXciB5Pk/4LwE3W9k3ARrfydSISLyKzgHnAh1YTUJeInCnOae2+4PaakBFJnbjurlyaR1xMFP/+Sg2vVTXR1t3vl+O0dPXx5d+Wc++GSlbMTOe1O8/liqWRu5iFmrymp0/hRzcsY3dTF999oSrY4XzMeIZsPomz0zZLROqA7wA/ANaLyJeAo8B1AMaYKhFZD1QDg8Ct1sgdgFv4aMjmq4RkJ24n09MSmJocH+xQAiotMZbbL5zLf/1xP389cByA+bnJnDFrKqtmZXLGrExyvFxw4vWqJu5+roKTfYPct7aYvz+r0Gdr0yoVii4oyuGW8+fwyNsHOGN2Jlcvzx/7RQEgIThU/mPKyspMeXl5QI51wYNvMy8nmUe/UBaQ44WavkE7FXUdbD7UxuZDbWw93EZ3v/Mze1ZWEqsKM1k1y/koyBzfoiEn+wb5/ovVPF1eS7EtlYfWLWPeJJnHX6mxDNod/M3/bKaivoMXv3Y2c3MCd+6LyFZjzCeSmSZ9S2fvACXffZ1/vmQ+X7tont+PFw4G7Q6qGzvZfND5IbDlcBsdPQMA5KVPOf0BcMasTGZlJX1iQYqtR9q48+md1J44xS3nzeHrF88nLkZvAleRpamjlzUPv8fU5Diev/XsgK2hPFLS1wnXLNUNzumUF0fQyJ2xxERHWQuDp/OVc2fjcBj2Nnex+WAbHx5q4719LWzY7hyWlpUczxnWh8DKwkxermjgkbcPMD19Cuv/8VOsLMwM8rtRKjimpSXw4xuWcdOvP+S+jVU8eN3SoMajSd8SqZ24ExEVJSyYlsqCaancdFYhxhgOtnbz4aE2Nh88zuZDbbxc0Xh6/+tK87nv8mKdM0dFvHPnZ3PbBc5+szNmZXJdWcHYL/ITTfqWivoObGkJZEVYJ643RIQ52cnMyU7mxlUzMMZQd6KHDw+1kZuawDnzdNy9Ui5fv3g+5YdP8K8bKynJT6doWnD6trSB1VJR38Gi6XqV7w0RoSAzkWtK8zXhKzVEdJTw0I3LSI6P5atPbKW7bzAocWjSxznC5FBrtzbtKKX8KiclgYfXLeNQazf/8nwlwRhIo0kfqKrvwBhYkh9Zd+IqpQLvrLlZ3HHRfDZsrw/KkqWa9HFfE1ev9JVS/nfbhXM5Z24W33mhiprGzoAeW5M+zpE7uanx5KR4d9epUkqNR3SU8JN1y0ibEsutT2zjZADb9zXpA5UNndqer5QKqKzkeB6+cTmHj3dzz3MVAWvfj/ik3903yIGWk9q0o5QKuDNnT+WfP1PEizsbeGLz0YAcM+KTfnVjp7MTV5O+UioIbjlvDufNz+b/vlQdkOnNIz7pV9TpnbhKqeCJihJ+fMMyMhPjuPX32+jsHfDv8fxaexiorO8gOyXe66mDlVLKU5lJcfz0b5ZTd6KHu5/d5df2/YhP+hURtiauUio0lRVmctfqIl6paOLx94/47TgRnfRP9WsnrlIqdNz86dlcuCCHf3u5ml117X45RkQn/eqGThzaiauUChFRUcJ/XreU7OR4bv39ttPrV/j0GD6vMYzodMpKqVCTkRTHTz+/gtIZGfhjRdGInlq5or6TrOR4clN1OmWlVOhYMSODFTMy/FJ3xF/pL8lL/cQyf0opNVlFbNLv6bezr7lLm3aUUhHFq6QvIneKSJWIVIrIkyKSICKZIvKGiOyzfma47X+PiOwXkT0istr78D1X3ejsxF2kSV8pFUE8TvoikgfcDpQZYxYD0cA64G7gLWPMPOAt63dEpNh6fhFwKfAzEYn2LnzPaSeuUioSedu8EwNMEZEYIBFoAK4Efms9/1vgKmv7SuApY0yfMeYQsB9Y5eXxPVZR38HUpDhsaXonrlIqcnic9I0x9cCDwFGgEegwxrwO5BpjGq19GoEc6yV5gPsyMXVWWVBU1newOC9NO3GVUhHFm+adDJxX77OA6UCSiPztaC8ZpmzYCSZE5GYRKReR8paWFk9DHFHvgJ19zSe1aUcpFXG8ad65GDhkjGkxxgwAzwFnAcdExAZg/Wy29q8DCtxen4+zOegTjDGPGmPKjDFl2dnZXoQ4vJrGTuwOo9MvKKUijjdJ/yhwpogkirON5CKgBngBuMna5yZgo7X9ArBOROJFZBYwD/jQi+N77HQnbr4mfaVUZPH4jlxjzGYReQbYBgwC24FHgWRgvYh8CecHw3XW/lUish6otva/1Rhj9zJ+j1TUd5CZFMd07cRVSkUYr6ZhMMZ8B/jOkOI+nFf9w+1/P3C/N8f0hYr6ThZN1ztxlVKRJ+LuyO0dsLPvmN6Jq5SKTBGX9Hc3dTHoMJr0lVIRKeKSfoXViasjd5RSkSjikn5lXQfpibHkZ0wJdihKKRVwkZf0G5xr4monrlIqEkVU0u8btLP3WJc27SilIlZEJf09TV0M2LUTVykVuSIq6Z/uxJ2uSV8pFZkiKulX1neQNiWWgkztxFVKRaaISvoV9R0s1jVxlVIRLGKSft+gnT1N2omrlIpsEZP09zad1E5cpVTEi5ikX9mga+IqpVTEJP2K+g5SE2KYkZkY7FCUUipoIibp65q4SikVIUm/f9DB7kbtxFVKqYhI+nuPddFvd2jSV0pFvIhI+qfXxNWkr5SKcBGR9CvqO0iJj2GmduIqpSJcRCT9yvoOFuWlEhWlnbhKqcg26ZP+gN1BTZOuiauUUhABSX/fsZP0D2onrlJKgZdJX0TSReQZEdktIjUi8ikRyRSRN0Rkn/Uzw23/e0Rkv4jsEZHV3oc/Nu3EVUqpj3h7pf8QsMkYswBYCtQAdwNvGWPmAW9ZvyMixcA6YBFwKfAzEYn28vhjqqjvIDk+hsKpSf4+lFJKhTyPk76IpALnAr8CMMb0G2PagSuB31q7/Ra4ytq+EnjKGNNnjDkE7AdWeXr88aqo76B4unbiKqUUeHelPxtoAX4tIttF5JcikgTkGmMaAayfOdb+eUCt2+vrrLJPEJGbRaRcRMpbWlo8DnDQ7qCmsVObdpRSyuJN0o8BVgCPGGOWA91YTTkjGO5S2wy3ozHmUWNMmTGmLDs72+MA9zWfpG/QoUlfKaUs3iT9OqDOGLPZ+v0ZnB8Cx0TEBmD9bHbbv8Dt9flAgxfHH9PpNXE16SulFOBF0jfGNAG1IlJkFV0EVAMvADdZZTcBG63tF4B1IhIvIrOAecCHnh5/PCrrO0iKi2Z2lnbiKqUUOJtovPE14AkRiQMOAv+A84NkvYh8CTgKXAdgjKkSkfU4PxgGgVuNMXYvjz+qyvoOFk1P005cpZSyeJX0jTE7gLJhnrpohP3vB+735pjjNWh3UN3Yyd+smhmIwymlVFiYtHfkHmjppnfAweK81GCHopRSIWPSJv0KvRNXKaU+YdIm/cr6DhLjopmdnRzsUJRSKmRM2qRfUd9BsS2VaO3EVUqp07wdvROylhWkY0tLCHYYSikVUiZt0v/XtcXBDkEppULOpG3eUUop9Uma9JVSKoJo0ldKqQiiSV8ppSKIJn2llIogmvSVUiqCaNJXSqkIoklfKaUiiBgz7IqFIUNEWoAjHr48C2j1UShal9aldWld4VTXTGPMJ9abDfmk7w0RKTfGDDffv9aldWldWlfE1OVOm3eUUiqCaNJXSqkIMtmT/qNal9aldWldWtdHJnWbvlJKqY+b7Ff6Siml3GjSV0qpCDIpk76IpIvIMyKyW0RqRORTXtRVJCI73B6dIvJ1D+u6U0SqRKRSRJ4UEY+X9hKRO6x6qjyJR0QeE5FmEal0K8sUkTdEZJ/1M8OLuq6zYnOIyLiHnY1Q1wPW/+UuEdkgIule1PV9q54dIvK6iEz3tC63574pIkZEsryI67siUu92nl3mTVwi8jUR2WP9H/yHF3E97RbTYRHZ4UVdy0TkA6uuchFZ5UVdS0XkfRGpEJEXRSR1nHUViMifrLxQJSJ3WOUTPvdHqWvC5/4odXl07o/KGDPpHsBvgS9b23FAuo/qjQaacN70MNHX5gGHgCnW7+uBv/cwjsVAJZCIc/WzN4F5E6zjXGAFUOlW9h/A3db23cAPvahrIVAEvA2UeRnXZ4AYa/uHXsaV6rZ9O/BzT+uyyguA13DeQJjlRVzfBb7pwbkwXF0XWOdEvPV7jjfv0e35/wTu8yKu14HPWtuXAW97UdcW4Dxr+4vA98dZlw1YYW2nAHuBYk/O/VHqmvC5P0pdHp37oz0m3ZW+9Yl/LvArAGNMvzGm3UfVXwQcMMZ4eodwDDBFRGJwJuwGD+tZCHxgjDlljBkE3gGunkgFxph3gbYhxVfi/MDE+nmVp3UZY2qMMXsmEtModb1uvU+AD4B8L+rqdPs1CRjXSIYR/r0Afgx8a7z1jFHXhI1Q1y3AD4wxfdY+zd7GJSICXA886UVdBnBdkacxzvN/hLqKgHet7TeAa8ZZV6MxZpu13QXU4Lwgm/C5P1Jdnpz7o9Tl0bk/mkmX9IHZQAvwaxHZLiK/FJEkH9W9jnGe9EMZY+qBB4GjQCPQYYx53cM4KoFzRWSqiCTivGoq8LAud7nGmEYr3kYgxwd1+toXgVe9qUBE7heRWuDzwH1e1HMFUG+M2elNPG5us77GPzbeprURzAc+LSKbReQdEVnpg9g+DRwzxuzzoo6vAw9Y//YPAvd4UVclcIW1fR0enP8iUggsBzbj5bk/pC6vjFKX1+c+TM6kH4Pzq+AjxpjlQDfOr2teEZE4nCfZHzx8fQbOq4lZwHQgSUT+1pO6jDE1OL/qvQFsAnYCg6O+aBIQkXtxvs8nvKnHGHOvMabAquc2D2NJBO7Fiw+NIR4B5gDLcF4U/KcXdcUAGcCZwF3AeutK3Rs34uEFj5tbgDutf/s7sb6Ne+iLwK0ishVnc0j/RF4sIsnAs8DXh3z7m7BA1OWrcx8mZ9KvA+qMMa5PyWdwfgh467PANmPMMQ9ffzFwyBjTYowZAJ4DzvI0GGPMr4wxK4wx5+L86uvNFZjLMRGxAVg/x9UsEAgichOwFvi8sRo4feD3jLNZYBhzcH6A7xSRwzi/dm8TkWmeVGaMOWaMsRtjHMD/AOPq5BxBHfCccfoQcOCcvMsjVnPk54CnvYgJ4Cac5z04L548fo/GmN3GmM8YY0pxfhgdGO9rRSQWZ2J9whjjisejc3+EujwyUl2+PvcnXdI3xjQBtSJSZBVdBFT7oGpvr3SOAmeKSKJ11XURznY7j4hIjvVzBs4/SG+vwgBewPmHifVzow/q9JqIXAp8G7jCGHPKy7rmuf16BbDbk3qMMRXGmBxjTKExphBnol1hnX+exGVz+/VqnM0XnnoeuNCqdz7OwQzezNZ4MbDbGFPnRR3gbMM/z9q+EC8uVNzO/yjgX4Cfj/N1gvMbRo0x5kduT0343B+lrgkbqS5fnvunedsTHIoPnF+Ry4FdOP8AMrysLxE4DqR5Wc/3cCaZSuB3WKMrPKzrPZwfZjuBizx4/ZM4mxEGcCasLwFTgbdw/jG+BWR6UdfV1nYfcAx4zYu69gO1wA7rMd4RN8PV9az1778LeBFnZ5lHdQ15/jDjH70zXFy/AyqsuF4AbF7UFQf8r/U+twEXevMegd8A/+SD8+scYKt1zm4GSr2o6w6cI1z2Aj/Aml1gHHWdg7NDeZfb+XSZJ+f+KHVN+NwfpS6Pzv3RHjoNg1JKRZBJ17yjlFJqZJr0lVIqgmjSV0qpCKJJXymlIogmfaWUiiCa9JVSKoJo0ldKqQjy/wG6JD4tFwJ6tgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "temp_2 = new4_2.iloc[new4_2.shape[0]-1,9:26]\n",
    "\n",
    "temp_2 = temp_2/dust_df_GJ['미세먼지 보정(㎍/㎥)'].mean()\n",
    "\n",
    "plt.plot(temp_2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "id": "6188fede",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "19.669715979254136"
      ]
     },
     "execution_count": 60,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dust_df_GJ['미세먼지 보정(㎍/㎥)'].mean()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 61,
   "id": "e56d9341",
   "metadata": {},
   "outputs": [],
   "source": [
    "new4_3 = bususer_GJ.copy()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 62,
   "id": "3a281844",
   "metadata": {},
   "outputs": [],
   "source": [
    "new4_3.loc['시간 별 평균 이용자 수',:]= new4_3.loc[:,'0':'23'].mean(axis=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 66,
   "id": "dfec7953",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x24b7808dfa0>]"
      ]
     },
     "execution_count": 66,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3dd3yUVfb48c9JQgi99yKIFAElQkSUYEMFdVUsFCsoLoplldXd1VX3u67r/lx1dVcsu1YQkWJD7AoiAtISTCBUQVroVVogJDm/P+4TGULqzGRmkjnv12temdyZuc8deHLmzrn3uVdUFWOMMdEhJtwNMMYYEzoW9I0xJopY0DfGmChiQd8YY6KIBX1jjIkiceFuQEkaNmyobdq0CXczjDGmQklNTd2pqo0Klkd80G/Tpg0pKSnhboYxxlQoIrK+sHJL7xhjTBSxoG+MMVHEgr4xxkQRC/rGGBNFLOgbY0wUsaBvjDFRxIK+McZEEQv6oTZnDqSmhrsVxpgoVWLQF5E3RWS7iGT4lCWKyDwRSRORFBHp6fPYwyKyWkRWikg/n/IeIrLEe+wFEZHgv50K4Lbb4M47w90KY0yUKk1PfwzQv0DZ08DjqpoI/MX7HRHpDAwBuniveVlEYr3XvAKMANp7t4J1Vn4HDsBPP8GPP7r7xhgTYiUGfVX9HthdsBio7d2vA2z27l8FTFTVI6q6FlgN9BSRZkBtVZ2rbquut4EBwXgDFcqSJaAKubmwYEG4W2OMiUL+5vTvB54RkY3As8DDXnkLYKPP8zK9shbe/YLlhRKREV7aKGXHjh1+NjECpaUduz97dvjaYYyJWv4G/ZHAKFVtBYwC3vDKC8vTazHlhVLVV1U1SVWTGjU6YZG4iistDerVg9NOs6BvjAkLf4P+UOBD7/57QP5AbibQyud5LXGpn0zvfsHy6JKeDomJkJwMc+dCTk64W2SMiTL+Bv3NwHne/QuBn7z7U4EhIlJVRNriBmwXqOoWYL+I9PJm7dwCfBxAuyue3FxYvPhY0D9wwOX4jTEmhEpcT19EJgDnAw1FJBP4P+C3wH9EJA44jJuVg6ouFZHJwDIgB7hbVXO9qkbiZgJVA77wbtHjp58gKwu6dYPevV3ZnDlwxhnhbZcxJqqUGPRV9foiHupRxPOfBJ4spDwF6Fqm1lUm+YO4iYnQujW0bOny+vfcE952GWOiil2RGyrp6VClCpx6Koi43v7s2W4KpzHGhIgF/VBJS4MuXSA+3v2enAybNsGGDeFtlzEmqljQD5W0NJfPz5ec7H7a1E1jTAhZ0A+Fbdtg61aXz8932mlQq5YFfWNMSFnQD4X0dPfTN+jHxsLZZ7sZPMZUcmlp8Prr4W6FAQv6oZE/c6dbNz77DGbM8MqTkyEjA/buDVvTjAmFu+6CESPcMJYJLwv6oZCW5qZp1qvHXXfBqFFeee/ebvbO3LlhbZ4x5Sk11Z3iqvDee+FujbGgHwppaZCYyObNbrLOkiXwyy/AWWe5NI/l9U0lNno01KgBnTrBpEnhbo2xoF/esrJg5UpITGT+fFeUlwfz5uH+Erp3t6BvKq3t22HCBPhf0qt8nnMJ8+Yp69eHu1XRzYJ+ecvIcFE+MZF589z1Wcd17nv3dmvrZ2eHtZnGlIfXXoPsbGXgz0/TdvU39GQBkyeHu1XRzYJ+efNZfmHePLfUTmKiz6Sd5GQ4fBgWLQpbE40pD0ePwiuvwL1J84jfuAaA3zWZbCmeMLOgX97S0qB2bXJatmHhQujVy3Xu581zfxS/Lr5mKR5TyUyZ4mbrjGo0DqpVg/PP58rDk1mUmsfq1eFuXfSyoF/evCtxl2QIWVku6Ccnu1R/WhrQtCm0a2fz9U2lM3o0tD8pmzbzJ8GAATB8OLV+yaQX8yzFE0YW9MtTXt6va+jPm+eK8nv64NO5T062xddMpZKWBrNmwdPnf47s3g033wxXXglVq3J/s0mW4gkjC/rl6eef3WYpXtBv3BjatIHmzaFt2wJBf+dOWLUqnK01JmhGj4bq1eGyXePciX/xxVC7Nlx6KZcdfI8li/NYsSLcrYxOFvTLk8+VuPPnu16+eLsFJye7jI4qx2+qYkwFt2sXvPsujBi4h/ivP4Xrr4c4b+uOQYOouW8Lycyx3n6YWNAvT2lpEBvL7mZdWLnSBf18ycluHbY1a3BXrTRoYIO5plJ4/XU3Ie2Bk95zU5FvvvnYg1dcAQkJ3N/cpXgsoxl6FvTLU3o6nHoqCxYnAMcH/ePy+vmbqlhP31RwOTnw8stw/vnQ8ttxbtOg7t2PPaFmTbj8cvrvf5+Vy3PJyAhbU6OWBf3y5C2/MG8exMRAUtKxh049FerV84nzvXu7nP727WFpqjHB8MknbqmRPw1a63o0N998LKeZb/Bgqu/fxvnyvaV4wsCCfnnZuRMyM6FbN+bNg65d3fL5+WJi4JxzCgzmgvX2TYX2wgtubcGLt7/jCm688cQnXXYZVK/O/c0nW4onDCzolxdvDf28092aO2eddeJTkpNhxQr3+UCPHlC1qgV9U2EtWQLffQd3jVRix4+D885znwAF1agBV1zBRfs+YO3qHH78MeRNjWoW9MuLF/TX1OzG3r3H5/Pz5Xfuf/gBF/DPPNMGc02F9eKLkJAAd3ZfAD/9dPwAbkGDBlFt/w76xnxnKZ4QKzHoi8ibIrJdRDIKlN8rIitFZKmIPO1T/rCIrPYe6+dT3kNElniPvSBSMNFXyaSlQYsWzFnVCCg86CcluX3Sj0vxpKbCoUOha6cxQbBnD7zzDtxwA9T55B0X/a+7rugXXHop1KzJ/c0nMXmypXhCqTQ9/TFAf98CEbkAuAo4XVW7AM965Z2BIUAX7zUvi0is97JXgBFAe+92XJ2Vjrf8wrx57pqUTp1OfEpCggv8xw3m5uTAwoUhbaoxgXrzTddXuffOozBxorv6tk6dol9QrRpceSUX7PmQzHVHWbAgdG2NdiUGfVX9HthdoHgk8JSqHvGekz/l5CpgoqoeUdW1wGqgp4g0A2qr6lxVVeBtYECw3kTEOXIEli//debOWWe5gdvCJCe7GJ+VhRvZBUvxmAolNxdeegn69IHErV+6QariUjv5Bg0i4eBu+sdNt7V4QsjfnH4HoI+IzBeRmSJyplfeAtjo87xMr6yFd79geaFEZISIpIhIyo4dO/xsYhgtWwY5OWR1SmTJksJTO/l693arbaakAPXrQ5cuNphrKpTPPoO1a+Hee4Fx46BhQ+jXr8TX0a8f1K7Nfc0mM3myW6rKlD9/g34cUA/oBfwBmOzl6AvL02sx5YVS1VdVNUlVkxo1auRnE8PIW35hSUwieXnFB/38zv1xKZ4ffnDdJ2MqgNGjoUULGHD+Xpg6FYYMcbsFlSQhAa66ij67PmJbZrZtFR0i/gb9TOBDdRYAeUBDr7yVz/NaApu98paFlFdOaWlQowYzNrQDCp+uma9hQ3eh1nGDub/8AkuXln87jQnQ8uUwbRqMHAlVPn7fpTZLk9rJN3gwVQ/t5fIq39gsnhDxN+hPAS4EEJEOQDywE5gKDBGRqiLSFjdgu0BVtwD7RaSX943gFuDjgFsfqdLS4PTTmbcghvbt3bI6xcnv3OflYRdpmQrlxRfdDLTf/haX2unQwU09Lq2LL4a6dflds8m89559wQ2F0kzZnADMBTqKSKaIDAfeBE72pnFOBIZ6vf6lwGRgGfAlcLeq5v83jgRexw3urgG+CPq7iQSqkJ6OdnODuMWldvIlJ7spb8uX49ZebtbMBnNNxPvlFxg71i2i2ThrPXz/feHLLhQnPh6uvpreO6awZ+thZs0qv/YaJ66kJ6jq9UU8dFMRz38SeLKQ8hSga5laVxGtXw+//MLu1ols3Vq6oO+7+FqXLnJsUxVjIthbb8HBg94A7vjxrvCmQsNC8QYNIv6tt7gi/msmTbqS888PZitNQXZFbrB5g7iLcrsBxefz87VrB02aFNgsfcMG2Lix2NcZEy55eW6a5tlnQ4/u6lI7ffq4b6pl1bcv1K/PfU0n8cEH7lIVU34s6AdbWhrExPDN1tNISIDTTy/5JVKwc2+bqpgI9+WXsHq118tPTXWLSPnTywc30+eaazhr+1T278hixoygNtUUYEE/2NLSoEMHZqVWJympdDPXwMX5tWth82agWze3KJUFfROhRo+Gpk3h2mtxvfz4eBg40P8KBw2iyuEDXJPwhc3iKWcW9IMtPZ3c0xJZtKh0+fx8x03aiYtz35str28i0KpVrqc/ciTEy1GYMMHtiFWvnv+VXnABNGzIvU0m8+GHbsMtUz4s6AfT3r2wbh2bGiWSnV22oJ+Y6DaSPi7Fs3gx7NtXLk01xl8vvui+wY4YAXz9NezYUba5+YWJi4NrryVp6ycc2XOQadOC0lRTCAv6weQtp/xjnhvELUvQr1LFDfoeN5iblwfz5gW5kcb4b/9+GDMGBg1y6R3eecddiHLppYFXPngwcUcOMbD655biKUcW9IPJm7nz5dZEWrZ0l6aXRXIy/Pij+8PirLMgNtZSPCaijB3rzs9778V9C50yBQYPdjn9QJ17LjRpwt2NJzNlittc3QSfBf1gSk+HJk34Kr1pmXr5+Xr3dp37+fNxeyt262ZB30SMvDyX2jnzTG8q8gcfuMjs76ydgmJj4brr6L7lM3L3HeCrr4JTrTmeBf1gSkvjSOdE1q4tW2on39lnuyWYj0vxzJ/vluE0JsymTYOVK+F3v/MKxo2DU07x72QvyqBBxB7J4vqan1qKp5xY0A+W7GxYupQN9cqez89Xu7ab13/cYO6hQ7+mjYwJp9GjoXFjb2bmxo1uQ9ybbirbsgslSU6GZs24q+Ekpk61TeTKgwX9YFmxArKzSc1JJC4Ounf3r5revWHuXO+qRLtIy0SINWvcuvl33OG2c2b8eLfOVLBSO/liYmDgQE7f9AUxB/fx+efBrd5Y0A8eb+bOl1sT6dbN7Qbnj+Rkt57J4sW4keC2bS2vb8LupZdcyv3OO3HBftw4txlEu3bBP9jgwcQePcJNtadaiqccWNAPlrQ0tFo1PszoEFCKM/8ireNSPLNn287RJmwOHHB74F57LTRvjks3LlsW+Nz8ovTqBS1bcmf9yXz2mTu+CR4L+sGSlkbWyV3Zfyg2oKDfsiW0bl1gMHfbNvj556A005iyeucdt4zycQO4VaoEtuxCcbwUT9fML4nP2ssnn5TPYaKVBf1gUIW0NNbXSwQCn8yQv/iaKoV0/Y0JHVU3gNu9u5tdRk4OvPsuXH55ybsDBWLwYGJyjjKs7seW4gkyC/rBkJkJu3eTkpNIgwaBpzmTk93Ca+vW4fZSrFvXgr4Ji2+/dZmce+/1JulMm+a+eZZXaidfz55w0kn8tu5kvvjCfdMwwWFBPxh8BnF79Qp8Bttxk3ZiYlyBzeAxYTB6tNvHecgQr2DcOLew2uWXl++BRWDQIE7d+DU1snfzceXdXDXkLOgHgzePfuq604JynUqXLlCnToHN0pcvh507A6/cmFJatw4++cTtf5uQgFt/4aOP3MI7VauWfwMGDSImN4fh9acweXL5Hy5aWNAPhrQ0DjU/hQPUCkrQj40tsLJyftf/hx8Cr9yYUsjfBlHELaEMuICflVX+qZ18PXrAySdze+1JfP2120faBM6CfjCkpbGuXiIibl2SYEhOhqVLvRP9zDPdglaW4jEh8PPPrtPx+efw3HPQqpX3wLhx7rqRc84JTUO8FE+HjdOpfXQnH30UmsNWdhb0A7V/P6xZQ8rRRDp3dmmZYMiftPPDD7jv1klJNphryt20aa6PkZkJX3zhM01z0yaYPj34yy6UZPBgJDeXOxp+aLN4gsSCfqAWLwbgi83dSrUJemmdeabbV+LXzn3v3pCSYuvNmnKhCs8/D/36QbNmsHAhXHKJzxPefdc9KVSpnXzdukH79txaczLTp7v9WkxgSgz6IvKmiGwXkYxCHntQRFREGvqUPSwiq0VkpYj08ynvISJLvMdeEAlld6EceYO4sw4kBnWxwerVXUrzuMHc7GwX+I0JoqwsuOUW+P3vYcAAt/bTCdOOx41z6ym3bx/axonA4MG02zCDBrnb+PDD0B6+MipNT38M0L9goYi0Ai4GNviUdQaGAF2817wsIrHew68AI4D23u2EOiuktDQO12zAJloENeiDi/MLFsCRIxzLo1qKxwTRhg3uPBs/Hp54At57z23lcJz0dFiyJPS9/HyDBiF5eYxsbCmeYCgx6Kvq98DuQh56Hvgj4LsozFXARFU9oqprgdVATxFpBtRW1bmqqsDbwICAWx8J0tNZVzeRmjWFzp2DW3Xv3i7gL1qEmyzdqZMN5pqg+f57N1S0ejV8/DE8+qi7LOQE48a5XOPgwSFvIwBdu0KnTgytNomZM2Hr1vA0o7LwK6cvIlcCm1Q1vcBDLYCNPr9nemUtvPsFy4uqf4SIpIhIyo5ITuLl5MCSJSzMTqRnTzfVMpjyZ2oel+KZM8dtYWSMn1Th5Zehb1+oX9/t03PFFUU8OTfX5fMvvdR1PMLBS/G02fA9jfO28P774WlGZVHmoC8i1YFHgL8U9nAhZVpMeaFU9VVVTVLVpEaNGpW1iaGzahUcPsy0Hd2CntoBt2FFhw4Fgv6ePe5CLWP8cOSIu9jq7rvdoO38+e4LZJG+/Ra2bAlfaiffoEGIKvc2+8BSPAHyp6ffDmgLpIvIOqAlsEhEmuJ68K18ntsS2OyVtyykvGLzBnEXaXAHcX3lr8Cgim2qYgKyeTOcfz688YZL5UydWoopxuPGuScV+VUgRDp3hq5duSl+ErNnuymlxj9lDvqqukRVG6tqG1Vtgwvo3VV1KzAVGCIiVUWkLW7AdoGqbgH2i0gvb9bOLUDFX00jPZ2c2HhW0Cmo0zV9JSfDrl1ub1LatYMmTWww15TZvHkuf79kCbz/vhu0LTR/7+vgQfjwQ7eEckJCSNpZrEGDaL1+Ns3ZxHvvhbsxFVdppmxOAOYCHUUkU0SGF/VcVV0KTAaWAV8Cd6tqrvfwSOB13ODuGuCLANsefmlprK/VldYnV6Fx4/I5xHGdexFbfM2U2RtvwHnnud3c5s1zm6GUaPdueOYZF/jDndrJN2gQAKNavGcpnkCoakTfevTooREpL0+1USOdUP1WveGG8j1Mw4aqw4Z5Bc89pwqqmzaV30FNpZCdrXr33e50ufhi1V27inlybq7qggWqjz+uevbZqjEx7oVJSe6xSNGtm2a2PltB9eefw92YyAakaCEx1a7I9dfWrbBjBz8cKr98PrjOff6mKsCx9Rmst2+KsX07XHSR29v2D39w6+jUr1/gSdu2uZz9jTe6tGHPnvDXv7oZO48+6tYAmTu3FHmgEBo8mBYb5tKKDbbypp8i6H+zgvEGcdNILLd8fr7evd1c6m3bgMREd7muBX1ThNRUl79fuNDNtnz6aTfNnqNHYdYseOQRd7l306buUtxp09yUzPHj3Uk2fz48/rhbdS0uLtxv53heiufB1pbi8VeE/Y9WIN7GKSviu5GYWL6H8u3cX3NNFXc5vA3mmkJMmQLXX++m+86ZA2c03AivfQlffumC+7597oKSc86BJ5+E/v1dRyKSevPFadcOevRg0I5J3PfjA6xa5aY1m9KrIP/TESgtjc1V23BKjzrEx5fvobp3d5MnjkvxpKXBgQPle2BToSxdCrfckMPwk6ax/PIHOePmrtC6NYwY4dbzGDTITd3ZudNdjvvnP7uTq6IE/HxDhtB0w0I6stJ6+36oYP/bkUN/TCPlaPnm8/PFx7t063Erbubmuq/hxuBW+L72GmVi3kBeXHkx1d8Y7ZbLfPZZyMhwi+y89pqbulO3bribG5gbb4SYGB5p9TYTJ4a7MRWPBX1/HDwIP60iNS80QR9c537RIjh0CJdrjYmxFI8B3IV7w4fDRT+9wmVHprjB2F274Jtv4IEH3P6blWRRW8B9mPXrx9X732bFslwyTlj/1xTHgr4/MjIQVdIIXdDv3dst9bNgAVC7Npx+ug3mGgD+8x9Y9l4Gz8c+4AZk//IXqFkz3M0qX8OGUXNvJn1lhvX2y8iCvj+8mTtbGnU7tpVcOTv7bNdZO27f3Llz3SeBiVqzZ8NjD2bxWe3riatfG956q3L16oty5ZVQty5/bDKGSZO8ZUpMqVjQ90daGvti6tCy90kh+/uqV8+tMHvcYO6BA7/u3GWiz7Ztbmz2pZp/4qR9GciYMW6+fTRISIAhQzh/94dsW73PLT9uSsWCvh+OpqSxKC+RXmeHtkeV37nPzcUu0opyOTluamavXZ9xyy+j4b77XGonmgwbRlx2FkNibc5+WVjQL6vcXCRjSUjz+fmSk90064wMoGVLNx1v5szQNsJEhMceg2UztvJu1WFufOepp8LdpNDr2RM6duS+2pbiKQsL+mW1Zg1xhw+yRLrRo0doD53fuf81xXPZZe76+l9+CW1DTFhNnQr/fCqP6S2HknD0AEyYEBmrYIaaCAwbRpc9s6myYTXz5oW7QRWDBf2y8gZxD7ZPpEaN0B66dWto0cIno3PbbW5X6wkTQtsQEzZr1riVE55t+R+6ZH4Nzz9P0PfprEhuugkV4bZYm7NfWhb0yyhvURpHiaPReaH/Qzth8bWkJDjtNLd2rqn0srLctVXd8n5k1LY/wVVXwR13hLtZ4dWyJXLxxfw2fizvTcpz412mWBb0y+jgD2ksozNJvauG5fjJybBxo7vAEhF3VU5Kis3iqeRU4a67YHX6QT6rewPSqBG8/np0TM8sybBhNMraQMdtM5k1K9yNiXwW9MsoZkk66ZTPnrilccKOiTfd5NZpsN5+pfbGGzBmDMzo/ntqZq6Et98O30blkWbAALR2bW6PG2MpnlKwoF8W27dTY+9mVlVLpH378DThtNOgVi2fFE+DBjBgALzzjtv12lQ6ixbBPffAX7t9xJmLXnUL5PftG+5mRY5q1ZDBg7mWD/jy/QMcPRruBkU2C/pl4S2nnNM1MWwLE8bFQa9eBabnDx/utrebMiU8jTLlZvdul8c/vX4mj2243a2D/8QT4W5W5Bk6lIScg5y/632+/TbcjYlsFvTL4PB8N3On/gXdwtqO5GSXwv91puZFF7mpPZbiqVTy8txMnS2ZuXzT/BZiso+4XVHKey3viuicc8hrdwrDY8fYhVolsKBfBnu/S2cDreh2YYOwtiM52Q3szZ3rFcTEwK23uk0y1q8Pa9tM8Dz1FHz2Gcy47BnqpM6A0aNtx5CiiBBz6zD65M4k9f21lukshgX9MohdkkY63ejZM7ztOOsst/nRcSmeW291P996KyxtMsE1bZq76vbRSxbQ6/PH3CI7w4aFu1mR7eabUREG7H+br78Od2MilwX90srKov6OFWQ2SKRevfA2pUYNOOOMAsvpn3SSS/O89RY2Wbliy8x06+p0b7+fx1ffgDRvDv/9r03PLEnr1ugFF3JrzFgmTcgLd2siVolBX0TeFJHtIpLhU/aMiKwQkcUi8pGI1PV57GERWS0iK0Wkn095DxFZ4j32gkjFOoM1YymxmkvuaeW8IW4pJSe7jbOOm6kwfLibwD99etjaZQKTne069YcPw7RT7yVm3Vo3MyvcPY0KIubWYbTJW8vOKbPdhkPmBKXp6Y8B+hco+wboqqqnA6uAhwFEpDMwBOjiveZlEYn1XvMKMAJo790K1hnRdk53M3fqXxgZQb93b3eF5o8/+hQOGAD169uAbgX2hz+4sZpvhk+kzpSx8Oij0KdPuJtVcVx9NTnVajIoawyffx7uxkSmEoO+qn4P7C5Q9rWq5u/eMQ9o6d2/CpioqkdUdS2wGugpIs2A2qo6V1UVeBsYEKw3EQq7Z6Sxn5p0vrxtuJsCHLtI67gUT9Wq7mKtKVPcdnmmQpk4EV54AR6/dR29xtzpds557LFwN6tiqVGDmCGDGCTvMWX8wXC3JiIFI6d/G/CFd78FsNHnsUyvrIV3v2B5oURkhIikiEjKjh07gtDEwMUtSWNJTDe6nh4ZwyDNmkG7dvDddwUeGD7c5QjeeScczTJ+WrYMbr8dzj0nh0dX3uTma44f7y7MMGUSM2woNfUA8Z99yP794W5N5AkogonII0AOMD6/qJCnaTHlhVLVV1U1SVWTGjVqFEgTgyMvjybb0tnWNDGi/gavucatrLxhg0/h6ae7hdjeeMMWGK8g9u93F2DVqAGf9HqSmB/muIHbtpHxrbLCSU4mq/nJ3HB0LFOnhrsxkcfvoC8iQ4HfADd6KRtwPXjfXWNbApu98paFlFcIR1auo2befvJOj4x8fr577nE/R48u8MDw4bBkiVuIzUS07Gy44QZYtQq+eHQOtf/9N7j5Zldo/BMTQ9URQ7mQb/n2LbtupSC/gr6I9Af+BFypqr5j5FOBISJSVUTa4gZsF6jqFmC/iPTyZu3cAnwcYNtDZt0UdyVugwvDeyVuQa1bux7ia6+57XJ/df31UK2aDehGuKNHYcgQ+PRT+N8/99L9XzdCmzbw4ovhblqFFzP0FmJQWs4Yx5494W5NZCnNlM0JwFygo4hkishw4EWgFvCNiKSJyH8BVHUpMBlYBnwJ3K2q+ZPGRwKv4wZ313BsHCDi7ZmRRi4xdLy2a7ibcoJRo9xyDMddk1WnDlx3ndtcxeatRaScHDfm/tFH8O/nldtTR7oJ+u++C7Vrh7t5FV+bNuzrcT435Y1lykeW5jyOqkb0rUePHhpuC5tfoSurdA53M4rUq5dqu3aqOTk+hd99pwqqY8eGrV2mcDk5qjfc4P57nnk6T/U//3G/PPlkuJtWqeS9+ZYq6P1nzg53U8ICSNFCYmpkTEWJcM23/cj2ZpGVz/c1apTbRu/TT30Kzz0XTjnFUjwRJi/PDbm8+y48+9gvPLjoBrjvPrff8Z/+FO7mVSoy8DqOVKlB55SxRMgkwIhgQb8EO5dtp3luJkdPD/Eu6GVwzTUuv//88z6FIm4P3e+/h59+ClvbzDF5eTBiBIwdC68Pn8sD4xLhvffgySfdbuexsSVXYkqvZk0OXnodg3QSH0+wNGc+C/olWPt+KgB1+0Zu0I+Lg3vvhZkzC1yhO3SoW4HzzTfD1jbjqJNHX6UAAB5eSURBVMLdd8Nbb+Ty9XlPctuYPu6DefZs+POfLeCXk3r3DaUO+9j6X9trIp8F/RIcmOmCfrtrIze9A+7Cnho1CvT2mzd3aYOxY93IoQkLVZfB+eS/maxqdREXz3wUGTTIfUKHa9/NKCHnn8feOidx5vKxbK4wk8TLlwX9EiQsS2VdlfbUblUn3E0pVt26LpszcSJs2eLzwPDhruCLCjNZqlJRhQcfhA2jp7AyoRsn717oNrsdP97NsjLlKyaGnBtu4WK+4cvXM0t+fhSwoF+C1ttT2dI8clM7vu67z3XoX3rJp/Dyy6FJExvQDQNVeOzBLE55biRTuJrqXdoiixa5tFvFWmS2Qmv4wFBiUI68YUuTgAX9Ym1dsoMWeRvJSawYQb9dO7jySncFf1aWV1ilittz79NPYevWsLYv2rx4ZwZDnjuTkfwXfeBB5IcfbOercGjXjo1t+nDBhjGsX2dz9i3oF+PnDxYBUC+CB3ELGjXKLbA5bpxP4W23uY1V3n47bO2KKqp8fvlL/PbVJFpV20neF18hzz5je9uGUcIdQ+nESmY9Oz/cTQk7C/rFOOgN4p58Xfcwt6T0zj3X7ar173/7rLfWqZNbi9kWYSt/O3eyqssALvv8Hla26EutnxcT0/+ScLcq6jW6ayCHY6qRMGlsuJsSdhb0i1FteSrr40+herOKM+Am4nr7y5fDV1/5PDB8uFvV67iNdaPIxo3u36BPH7j/fjeQumqVmzwfLDNmsL9dN05a/iVvd/83Xdd9SkzTxsGr3/ivdm3WJl5D350TWZ1xONytCSsL+kVQhdY7UtnWouKkdvINHuzW2z9u+ubAgVCzZvQN6O7bB4884nLp48e7ke7XXnML33TsCA0awMUXu7nyH30EmzaV/RhHj8Ijj6B9+7JpXy0eu3g+N8y/j9g4G6yNJA0fHEY99pL+RJSvt1zY2gyRdAvX2jvrF+1UBf1hwD/DcvxA/f3vbjmXjAyfwttvV61eXfWXX8LWrpDJzlZ9+WXVRo3cP8SNN6quW+ceO3pUNT1d9fXXVe+4Q7V7d9W4OPc8UG3WTPXKK90/4ldfqe7aVfRx1qxRPessVdBXuV2vv/KAZmeH5i2aMsrJ0a1VW+n3NS8Nd0tCgiLW3gl7UC/pFq6gP/ORr1RBV7w0LSzHD9SOHaoJCS7O/2ruXPdf/uqrYWtXucvLU506VbVTJ/dezztPdeHCkl936JD793nhBdWbblLt2PHYhwC4Fe2uv171uedUZ81SPXBAdfx41Vq19HC1Onodk/WKK1SPHCn3d2gCsPCSP2sOMbri203hbkq5s6BfRp+f+w9V0MNbdofl+MEwYoRq1aqq27d7BXl5qp07u55pZZSSonr++e607thR9eOP3Xv21969qtOnq/6//6d6zTWqrVod+xCIiVEF3XpKb23NOr3sMtXDh4P3Vkz52DFnpSrol32fDndTyp0F/TKa2eha3Vj15LAcO1iWLXP/w3/7m0/hv/6lJ+Z9Krj1613vHFQbNlR96SUttxzLli3um8Rjj+n8m17QOI7qJZeoZmWVz+FM8GXUOVtXxXfWvNwAOgQVQFFB3wZyC5GXB613LqqQg7i+Tj0V+vd3V+geOeIV3nyzu2CrMgzo/vILPPywG6R9/313f/VquOsu9x6D5MABN9Hnu+9g/PSmPLPiCu7a+TfOfvdezr0wjilTICEhaIcz5Wz3FcNon72MlRNSw92UsIigbb4jx88puzlF17LjjDvC3ZSAjRoF/fq5NXmGDgUaNXKX7Y4bB089VTEvGDp6FF59Ff76V9i5032Q/f3vbn3pMjhyxC1LtGkTbN58/M23bP/+E19bowZccYWbEFStWnDelgmNLo8PIuud+9j9rzFwY1K4mxNyFvQLsf6jRZwCNLikYvf0wc1G7NLFTd+85RZvyZfhw+GDD9wa7tddF+4mlp6qa/Mf/+i63uefD88+Cz1K/n86dMh945k+/Vgw37XrxOfFx7vFSZs3h9NOcx+YLVocK8u/2Y6GFVf9k+sys9kAuqW/ix7+F5JQNdxNCikL+oU4NMt97Ws9oOJciVsUEXct0m9/69ITF1wAXHIJtGzpUjwVJeinpMADD7hNYTp1gk8+cYvJlbBwWXY2vP46PPGEW3ro9NOhbVt3gXLz5icG9AYNbC20aHD0xmHUfXYiq57/lA4PXxvu5oRWYYn+SLqFYyB3esOBuqlqm5Aft7wcOuTGN6+4wqfw0UdVRVQ3bAhbu0pl/Xo3xx7cnPuXX3bz7EuQk6P69tuqbdu6lyYnq37/fQjaayqEvbtyNJPmmtH2N+FuSrnBBnJLJycH2uxKZXvLip/ayVetGowc6Rba/HXnxFtvdemSMWPC2bSiqbrNX7p2damoP//ZDdKOHOm2CivmZR9/DN26uXRWnTrw+efuC0KfPiFsv4loderHMr/9zXRc+wW5m7eFuzkhZUG/gFXz93Cy/ox2rzxBH45NaHnhBa/g5JPhwgvdVorBXH8mGHbudGmnYcMgMRGWLXP7yJaQSP/2Wzj7bBgwwKV1Jk6E1FS49FJL2ZgTJdw5lDhyWffkeP8rycuDw4fdch+7drn7Ea7EoC8ib4rIdhHJ8CmrLyLfiMhP3s96Po89LCKrRWSliPTzKe8hIku8x14Qicw/ww1T3HLKDfpVrqDftClcfz289Rbs3esVDh8O69bBjBnhbNrxvvjCjaB+8gn885+ubW3bFvuShQvdgHXfvm7WzWuvuc+JwYPdFsHGFOa8O05lYUxPGo55xs1o69/fnUR9+sBZZ7nlart2dVOC27Rxgz6NGrmvj9WquW+csbHufp060LChm0G2enW431qxSjOQOwZ4EfBdjP0hYLqqPiUiD3m//0lEOgNDgC5Ac2CaiHRQ1VzgFWAEMA/4HOgPRNwefllzXNBveWXlCvrgpm+OHeuC4h/+AFx9tdtn8Y033MkeTocOuUa9/LKbbvTlly5HU4xly+DRR906aQ0bwnPPueyPzZk3pVGjBnx3ziNUnfMY8dMziaseT5Ua8STUjCehTk2qN4ynSvV49xU5Pt7dirsfGwt/+5tb3HDu3Mg9EQtL9Be8AW2ADJ/fVwLNvPvNgJXe/YeBh32e9xVwtvecFT7l1wP/K82xQz2Q+3X9wbql6kkhPWYoXXCBW03g17HQu+92azXsDuNyEwsWqHbo4EZcf//7Ei9vXbtW9ZZb3EoItWqpPv54dKwhZ4Jv5UrVkSNV+/ZVbd1aj1tuKX/tvfPOc2tYPf206pQp7kr3IpfcmDrVvfDOO0P5NgpFEQO5/k7ZbKKqW7wPjS0ikr9oeAtcTz5fpld21LtfsLxQIjIC962A1mW84CYQ2dnQdncqO9r3oGnIjhpao0a5b7IffODSHwwf7iawjx8P99wT2sbk5MA//uF6R82auUn0F15Y5NO3bnWp/f/9z6VtRo2Chx5yvXxj/NGhg/tyme/QIVizxl0GsmqVm/iwapWbHLBjx7HnxcS4jE/79q6O/FviWVfQ+MEH3fUj557rcqqRprBPgoI3Tuzp7y3w+B7v50vATT7lbwDXAmcC03zK+wCflObYoezp//jdXlXQ9MFPhuyYoZabq9q+fYE11844QzUxMbQNWbXq1yWJ9YYbiv2msWeP6sMPu1WhY2PdQnIbN4awrcaoO0Xnz1d95x3Vv/xFdcgQtyp3zZrHvhlUq6b61afZquec4x5YsSJs7SXIUza3iUgzAO/ndq88E2jl87yWwGavvGUh5RFl48cun9+oX8W/KKsoMTFw330wf75LOwKut5+WBosWlX8DVF1XPTERVq6ECRPct4x69Qp96osvunHc//f/4Kqr3I5g//ufu7bMmFCqVw969oQbb4THH3enbmqqm7izZYubc9ChA1xxTRW+GT4RqlaFQYMgKyvcTT9eYZ8EBW+c2NN/BnjIu/8Q8LR3vwuQDlQF2gI/A7HeYwuBXoDgBnAvK82xQ9nTn9zzGVXQvG3bS35yBbZ/v2rduqoDB3oFu3e7vP6FF6p+9JFbjL88bN2q+pvfuC5R377Fdtf37FG9+mr31IsvVk1LK58mGRNMu3apJiW5PXlm/flzdwIft6lF6ODv0srABGALx/Lyw4EGwHTgJ+9nfZ/nPwKswQ32XupTngRkeI+9CEhJx9YQB/0v6g7RbQmtQna8cPrjH91AaP5mUvr3v7tdV/K/p3bponrXXaoTJ6pu3hz4AadMcZcFV62q+u9/uzxTERYudFfSxsW5laADWRLfmFDbu9dld2JjVTOueMj9PY0bF/J2+B30w30LVdA/dEh1Je11aYcBITleuG3Y4E7KBx7wKTx8WHX2bNV//EO1X7/jk5Xt26sOH646dqybPlNa+/a514EbN1i6tMin5uWpjh6tGh/vZhj98IPfb8+YsNq/3+3nE8dR3dKhjxuQWrYspG2woF+CBd+4Qdyl1z8RkuNFgiFDVGvXdnG5UEePum73s8+6PWPr1Tv2IdC6tdu45NVX3by3wrrjc+aonnyyW+Pn4YeL3Utw716XbgLVyy9X3bkzOO/RmHA5eNClJpuTqYdqNnLfng8cCNnxLeiX4IPffacKuuPtz0NyvEgwf747A/7zn1K+IDdXdfFi1x0fOFC1ceNjHwJNmriy0aNdAv7Pf3b5ozZt3J6yxUhNdVvQxsa6udDFZH6MqVCystww1sV8pXkiqsOGhezYFvRLMD7JbSOYt3VbSI4XKc45x3XGc3L8eHFenpuS9uqrbiXMli2PfQiA6q23FnvVVF6eWzQzPl61RQuXWTKmsjlyxG2x/DcedX8Xb70VkuMWFfTFPRa5kpKSNCUlpdyP82mdGzkr+3saZW0s92NFkvffd1eNf/SRW6gsIKqwfr1b0rJ5c7jooiKfum8fjBgBkya5JU/GjbOLrEzllZMDQ2/KZfiki+kTN4+4RQuQ07qW6zFFJFVVT9gazII+bju8TbU7QceOdFrxcbkeK9Lk5MApp8BJJ8HMmaE5Znq6+6BZs8btcvinP9nCaKbyy82F39+wlYcnJxLToB6N1i5EatUst+MVFfTtTw1YPGc/HViFJFW+RdZKEhcHv/ud65yX97VZqm6xt1693GbjM2a4vcwt4JtoEBsLz09oyoTfvEuDXatI7TkSzQt9p9v+3IDMT34kBqVx/+gL+uAuyK1Z02068re/uWB86FBwj3HggNu/fMQIt3JtWppbmsSYaBITA/dPvZCve/0fSSveYfyFb4R8OwsL+sCRuW5P3HoXRWfQr1MH/vtf1xP561/dmmd16rglxR98EKZMcfua+GvJEkhKcpetP/GEWzW5ceOSX2dMZSQC/Wc9wk8nXcS1M+/lr9csJjc3hMe3nD58XOsmknNm0CBrU7kepyLYuxd++AFmz4ZZs2DBArf6KLj9yPv0geRk97NNm+J3pFJ1m7bcc4/7EHn3XW9jdmMMum07+09JZMuBWjw9MIX/vVuruJ1Ay8wGcouwZw9sqd+ZKp1Oof3yqeV2nIrq8GG3qNSsWe6DYM6cYztvNW9+7AMgOdlteBUb6x47eNBt0fj22+6bw/jxbvcuY4yPmTPJu+BCJukgPrjmXd6dIMTHB6fqooJ+ED9XKqa02Qc4jxWsPXNwuJsSkRISoHdvdwO3JejSpce+CcyaBZMnu8dq14ZzznG3iRPdipj/93/w2GPHPgyMMT7OO4+Yvz/B9Y88wswPz+O66+5k8uTy3XQr6nv679w5m5v+14f9735Cret/U27HqczWr3cfAvkfBEuXupz9+PHFTtU3xoDrSV1+OTnTZnBmzlwaX3IGH30E1asHVq319IuQ7Q3i1jo/Ogdxg+Gkk9ztxhvd73v2uJ5KtWrhbZcxFUJMDIwbR1xiIjOzB9L661R+85s6TJ3qZtUF/XDBr7JiqbM6lT0JTd12fSYo6tWzgG9MmTRsCJMmUXv3OtLPvJ2Z3yn9+7sr14MtqoP+9u3Q6VAqu9taL98YE2a9e8M//sFJC99n4a0vkZNDuUzljOqg/+Psg3RiBbE9LegbYyLAgw/C5ZfTfdzv+eGFlMJ2EQ1YVAf9zZ+nEUsejS+1oG+MiQAxMTB2LDRtSsyQQcfmRwfzEEGvsQLJnucWm6mebEHfGBMhGjRw86DPOadcFqaK2tk7qlB3TSp7E5pQt3nzcDfHGGOO6dXL3cpB1Pb0N22CzodT2Xtyj+LXEjDGmEokaoP+j3MO0ZllxJ1lqR1jTPQIKOiLyCgRWSoiGSIyQUQSRKS+iHwjIj95P+v5PP9hEVktIitFpF/gzfff5i/SiSWPRv26h7MZxhgTUn4HfRFpAfwOSFLVrkAsMAR4CJiuqu2B6d7viEhn7/EuQH/gZREJ24osOfPdlbhVz7GevjEmegSa3okDqolIHFAd2AxcBYz1Hh8L5O+8ehUwUVWPqOpaYDXQM8Dj+0UV6v2cyr6ERtCyZTiaYIwxYeF30FfVTcCzwAZgC/CLqn4NNFHVLd5ztgD522W0AHx3Hc/0ykJu3Tromm2DuMaY6BNIeqcervfeFmgO1BCRm4p7SSFlhS7xKSIjRCRFRFJ27NjhbxOLtGhOFp1ZRpVeltoxxkSXQNI7FwFrVXWHqh4FPgTOAbaJSDMA7+d27/mZQCuf17fEpYNOoKqvqmqSqiY1atQogCYWbvOXi4kjl4ZRuieuMSZ6BRL0NwC9RKS6iAjQF1gOTAWGes8ZCnzs3Z8KDBGRqiLSFmgPLAjg+H7LXeAGcavYdE1jTJTx+4pcVZ0vIu8Di4Ac4EfgVaAmMFlEhuM+GAZ6z18qIpOBZd7z71bVEG4H7OTlQf11qexPaEitVq1KfoExxlQiUbdz1sqVkNUpkYZdmtAy46ug1WuMMZGkqJ2zou6K3NQ5h+nCUuJtENcYE4WiLuhv/XoxVcihwSUW9I0x0Sfqgn7uQjeIaxunGGOiUVQF/ZwcaLQ+lYMJ9d1O3sYYE2WiKugvXw6n5y7il1PsSlxjTHSKqqCf+sMRupJB1bMttWOMiU5RFfS3frOEeI5S7yIL+saY6BRVQT8vxQ3ixiTZGvrGmOgUNUH/yBFovCGVQwn1oG3bcDfHGGPCImqCfkYGnKGp7D+luw3iGmOiVtQE/dQfjnAaS2ynLGNMVPN7wbWKZtv0DOI5SpULLegbY6JX1PT081IWASBJFvSNMdErKoL+oUPQdHMqWQl14eSTw90cY4wJm6gI+unp0F1TOdDeBnGNMdEtKoJ+6txsTmcxCefY/HxjTHSLioHcbd8upSrZxJ9n+XxjTHSLip4+qe5KXBvENcZEu0of9PfvhxZbUzlctTa0axfu5hhjTFhV+qC/aBF0J5WDHbpDTKV/u8YYU6xKHwVT5x3ldBZTLdlSO8YYU+kHcrd+u4wEjkAfC/rGGBNQT19E6orI+yKyQkSWi8jZIlJfRL4RkZ+8n/V8nv+wiKwWkZUi0i/w5pcs5kc3iEsPC/rGGBNoeuc/wJeq2gnoBiwHHgKmq2p7YLr3OyLSGRgCdAH6Ay+LSGyAxy/W7t3QakcqR6rWglNOKc9DGWNMheB30BeR2sC5wBsAqpqtqnuBq4Cx3tPGAgO8+1cBE1X1iKquBVYDPf09fmmkpkIPUjnU4QwbxDXGGALr6Z8M7ADeEpEfReR1EakBNFHVLQDez8be81sAG31en+mVnUBERohIioik7Nixw+8GpszLoRvpNohrjDGeQIJ+HNAdeEVVzwAO4qVyilDYojda2BNV9VVVTVLVpEaNGvndwO3fLaMah0nobUHfGGMgsKCfCWSq6nzv9/dxHwLbRKQZgPdzu8/zW/m8viWwOYDjl8gGcY0x5nh+B31V3QpsFJGOXlFfYBkwFRjqlQ0FPvbuTwWGiEhVEWkLtAcW+Hv8kmzbBm33pJJdtSZ06FBehzHGmAol0Hn69wLjRSQe+Bm4FfdBMllEhgMbgIEAqrpURCbjPhhygLtVNTfA4xcpJQW6s4isjmcQb4O4xhgDBBj0VTUNSCrkob5FPP9J4MlAjllaqfNzeJA04pLvCMXhjDGmQqi0V+Rum7mC6mRBL1tD3xhj8lXavEe9n20Q1xhjCqq0Qf+JAalojRrQsWPJTzbGmChRadM7sigVEhMhtlxXejDGmAql0gZ9zjoLWrUq+XnGGBNFKm/Qf+65cLfAGGMiTqXN6RtjjDmRBX1jjIkiFvSNMSaKWNA3xpgoYkHfGGOiiAV9Y4yJIhb0jTEmiljQN8aYKCKqhe5YGDFEZAew3s+XNwR2BqkpVpfVZXVZXRWprpNU9YT9ZiM+6AdCRFJUtbD1/q0uq8vqsrqipi5flt4xxpgoYkHfGGOiSGUP+q9aXVaX1WV1WV3HVOqcvjHGmONV9p6+McYYHxb0jTEmilTKoC8idUXkfRFZISLLReTsAOrqKCJpPrd9InK/n3WNEpGlIpIhIhNEJCGAdt3n1bPUn/aIyJsisl1EMnzK6ovINyLyk/ezXgB1DfTalicipZ52VkRdz3j/l4tF5CMRqRtAXU949aSJyNci0tzfunwee1BEVEQaBtCuv4rIJp/z7LJA2iUi94rISu//4OkA2jXJp03rRCQtgLoSRWSeV1eKiPQMoK5uIjJXRJaIyCciUruUdbUSkRleXFgqIvd55WU+94upq8znfjF1+XXuF0tVK90NGAvc7t2PB+oGqd5YYCvuooeyvrYFsBao5v0+GRjmZzu6AhlAddzuZ9OA9mWs41ygO5DhU/Y08JB3/yHgnwHUdSrQEfgOSAqwXZcAcd79fwbYrto+938H/NffurzyVsBXuAsIGwbQrr8CD/pxLhRW1wXeOVHV+71xIO/R5/F/AX8JoF1fA5d69y8DvgugroXAed7924AnSllXM6C7d78WsAro7M+5X0xdZT73i6nLr3O/uFul6+l7n/jnAm8AqGq2qu4NUvV9gTWq6u8VwnFANRGJwwXszX7WcyowT1UPqWoOMBO4uiwVqOr3wO4CxVfhPjDxfg7wty5VXa6qK8vSpmLq+tp7nwDzgJYB1LXP59caQKlmMhTx7wXwPPDH0tZTQl1lVkRdI4GnVPWI95ztgbZLRAQYBEwIoC4F8nvkdSjl+V9EXR2B77373wDXlrKuLaq6yLu/H1iO65CV+dwvqi5/zv1i6vLr3C9OpQv6wMnADuAtEflRRF4XkRpBqnsIpTzpC1LVTcCzwAZgC/CLqn7tZzsygHNFpIGIVMf1moKxC3wTVd3itXcL0DgIdQbbbcAXgVQgIk+KyEbgRuAvAdRzJbBJVdMDaY+Pe7yv8W+WNrVWhA5AHxGZLyIzReTMILStD7BNVX8KoI77gWe8f/tngYcDqCsDuNK7PxA/zn8RaQOcAcwnwHO/QF0BKaaugM99qJxBPw73VfAVVT0DOIj7uhYQEYnHnWTv+fn6erjeRFugOVBDRG7ypy5VXY77qvcN8CWQDuQU+6JKQEQewb3P8YHUo6qPqGorr557/GxLdeARAvjQKOAVoB2QiOsU/CuAuuKAekAv4A/AZK+nHojr8bPD42MkMMr7tx+F923cT7cBd4tIKi4dkl2WF4tITeAD4P4C3/7KLBR1Bevch8oZ9DOBTFXN/5R8H/chEKhLgUWqus3P118ErFXVHap6FPgQOMffxqjqG6raXVXPxX31DaQHlm+biDQD8H6WKi0QCiIyFPgNcKN6Cc4geJdSpgUK0Q73AZ4uIutwX7sXiUhTfypT1W2qmquqecBrQKkGOYuQCXyozgIgD7d4l1+8dOQ1wKQA2gQwFHfeg+s8+f0eVXWFql6iqj1wH0ZrSvtaEamCC6zjVTW/PX6d+0XU5Zei6gr2uV/pgr6qbgU2ikhHr6gvsCwIVQfa09kA9BKR6l6vqy8ub+cXEWns/WyN+4MMtBcGMBX3h4n38+Mg1BkwEekP/Am4UlUPBVhXe59frwRW+FOPqi5R1caq2kZV2+ACbXfv/POnXc18fr0al77w1xTgQq/eDrjJDIGs1ngRsEJVMwOoA1wO/zzv/oUE0FHxOf9jgEeB/5bydYL7hrFcVZ/zeajM534xdZVZUXUF89z/VaAjwZF4w31FTgEW4/4A6gVYX3VgF1AnwHoexwWZDGAc3uwKP+uahfswSwf6+vH6Cbg0wlFcwBoONACm4/4YpwP1A6jrau/+EWAb8FUAda0GNgJp3q20M24Kq+sD799/MfAJbrDMr7oKPL6O0s/eKaxd44AlXrumAs0CqCseeMd7n4uACwN5j8AY4M4gnF/JQKp3zs4HegRQ1324GS6rgKfwVhcoRV3JuAHlxT7n02X+nPvF1FXmc7+Yuvw694u72TIMxhgTRSpdescYY0zRLOgbY0wUsaBvjDFRxIK+McZEEQv6xhgTRSzoG2NMFLGgb4wxUeT/A/5ZWaCQJW3rAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "temp_3 = new4_3.iloc[new4_2.shape[0]-1,9:26]\n",
    "\n",
    "plt.plot(temp_3, c = 'blue')\n",
    "plt.plot(temp_2, c = 'red')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "60bcbe38",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
